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ABSTRACT 


Sound radiation ultimately arising from gearboxes is a major 
contributor to helicopter and turboprop cabin noise. Application 
of the transfer function approach to predict the resulting 
interior noise contribution requires gearbox vibration sources 
and paths to be characterized in the frequency domain. In pre- 
vious work [J.Acoust. Soc. Am., 63 (1978), pp. 1409-1430], the 
frequency domain characterization of the vibratory excitation of 
a generic pair of meshing gears was shown to require the Fourier- 
series coefficients of the static transmission error of that 
pair. Static transmission error spectrum contributions at tooth- 
meshing harmonics were shown to be attributable to elastic tooth 
deformations and the average deviation of gear teeth faces from 
perfect involute surfaces with uniform spacing, whereas spectrum 
contributions at other gear rotational harmonics were shown to be 
attributable to deviations of individual tooth faces from the 
average face of uniformly spaced teeth. The present work con- 
tinues this analysis by representing tooth-face deviations from 
perfect involute surfaces in terms of Legendre polynomials which 
may be directly interpreted in terms of tooth-spacing errors, 
mean and random deviations associated with involute slope and 
fullness, lead mismatch and crowning, and analogous higher-order 
components. The contributions of these components to the spectrum 
of the static transmission error is discussed and illustrated 
using a set of measurements made on a pair of helicopter spur 
gears. The general methodology presented is applicable to both 
spur and helical gears. 
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INTRODUCTION 


For some time, it has been recognized that sound radiation 
ultimately arising from the gearbox is a major contributor to 
helicopter interior noise. Many analyses of interior noise mea- 
surements have shown that the dominant frequencies are those 
associated with the meshing action of the gears - i.e., the so- 
called tooth-meshing harmonics and their sidebands. Hence, one 
approach to reducing helicopter interior noise is to reduce the 
vibratory excitation caused by the imperfect meshing action of 
the gears. Such an approach requires a thorough understanding 
of this source of vibratory excitation - which is the subject 
of the present report. 

Consider a pair of -rotating gears that transmit torque, and 
assume for simplicity that one member of the pair rotates at an 
exactly uniform speed. The mating gear will have a small jitter 
superimposed on its mean rotational speed, where this jitter 
arises from the nonconjugate action of the imperfect loaded teeth. 
This unsteady motion component is the dominant ultimate source 
of the vibratory excitation and noise arising from the gear- 
box [i]. 

The static transmission error describes this source of 
vibratory excitation. In Ref. 1, a method is explained for 
writing the equations of motion of gear systems using the static 
transmission error as the input. Such equations are most con- 
veniently solved in the frequency domain. Since the static 
transmission error is periodic, it is described in the frequency 
domain by its Fourier-series coefficients, which therefore are 
of primary interest. 

Expressions for these Fourier-series coefficients are given 
in Ref. 1 in terms of the mean surface of the tooth faces on a 
gear, deviations of the individual tooth-face surfaces from this 
mean surface, and a general tooth-pair stiffness characteriza- 
tion. These Fourier-series coefficients - and the corresponding 
power spectra of the individual tooth— face deviation components — 
are expressed in terms of general series expansions of the tooth- 
face surfaces and the tooth-stiffness characterization. The 
present report extends the work carried out in Ref. 1 by further 
developing these series expansions in terms of general ortho- 
gonal functions, and then specializing the results to expansions 
in two-dimensional Legendre polynomials, which are uniquely 
suited for the description of the deviations of tooth faces from 
perfect involute (conjugate) surfaces. 
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After completion of these general results, the Fourier-series 
coefficients and power spectra of the various static transmission 
error components are developed from a set of tooth-spacing 
accuracy measurements and profile measurements made from a pair 
of helicopter spur gears. Some of these results are then com- 
pared with results that would have been obtained with helical 
gears of comparable quality — bringing into focus the very sub- 
stantial noise reduction achieved through the use of helical 
gears . 

Acknowledgements . I am indebted to Bell Helicopter Textron 
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measurements for the Applications portion of this work. The 
cooperation of Mr. Charles E. Braddock, Senior Transmission 
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Fischer of BBN. Discussions with Dr. Fred R. Kern have been 
helpful in various stages of the work. The report was typed by 
Ms. Carolyn Gianino and the art work was done by Ms. Laura 
Selvitella . 



HARMONIC ANALYSIS OF STATIC TRANSMISSION ERROR 
USING ORTHOGONAL EXPANSIONS 


Orthogonal Expansions in Tooth Coordinates 

Consider the expansion of an arbitrary function P^(y 5 z) in 
a set of orthogonal functions y^ m (y.,z); i.e., C 

V y ’ z) = . (i) 

m 


where y,z are the Cartesian tooth coordinates described in Sec. 
I.C of Ref. 1 and illustrated by the tooth shown in the upper 
right-hand corner of Pig. A of the same reference. In later 
sections, the expansion of Eq. (1) will be used to represent 
local tooth-pair stiffnesses K TC (y,z) and stiffness-weighted 
deviations of tooth faces from perfect involute surfaces. We 
assume that the expansion functions Tp„(y,z) have been 
normalized so that their mean— square values are unity — i.e.. 


1 

A' 



(y,z)dydz 


mm 


( 2 ) 


where A* denotes the area of the region in tjhe y,z plane where 
the expansion functions are defined, and is Kronecker's 

delta. Since A' has the dimension of length squared, we see from 
Eq. (2) that the functions T Cm (y,z) must be dimensionless. Thus, 
the expansion coefficients a m in Eq. (1) have the same dimension 
as F c (y,z) . 

We multi P!y Ec J* (1) b y Tom' (y,z)/A' , integrate both sides 
of the resulting expression over Q c , and use Eq. (2), we obtain 


a 


m 



(y,z)'F Cm (y,2)d y dz, 


( 3 ) 


which determines each expansion coefficient a m as a function of 
F c (y,z). Squaring Eq. (1), we obtain 
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Then, dividing Eq. (4) by A', integrating both sides of the 
resulting expression over and using Eq. (2), we have 


1 

A’ 


j Fq (y , z )dydz 



(5) 


provided that the expansion functions form a complete set. 
According to Eq. (5), for each value of m, a^ directly measures 
the contribution of the term a m 'i , Q m (y,z) in Eq. (1) to the mean- 
square value of Fc(y,z) defined by the left-hand side of Eq. (5). 
This property, which is a consequence of the normalization of 
Eq. (2), greatly facilitates interpretation of numerical values 
of the expansion coefficients. 


Expansions in Legendre polynomials 

We shall see later that the Legendre polynomials are a par- 
ticularly useful set of functions to use in our orthogonal expan- 
sions. For the expansion functions in Eq. (1), let us take 

' 1 'cm (y > z) = V (y), L* (z) ’ <6) 

where 


, ^ A 

yk y = 

(2k+l)^P k (2y/F), 

- (F/2 ) <y< (F/2 ) 

(7) 

( \ A 
2 k (z) = 

( 2 £+1 P^ ( 2z/D ) , 

-(D/2)<z<(D/2) , 

(8) 


where P^(x) and P£, ( x ) are Legendre polynomials as defined, for 
example, in Ref. 2. Thus, for the H'c m (y,z) defined by Eqs . (6) 
to (8), the expansion region fic is the rectangular domain y,z 
indicated in Eqs. (7) and (8). 

Equations (6) to (8) imply a mapping of the Legendre indices 
k and l onto the single index m. One such mapping is illustrated 
in Fig. 1. From the normalization of the Legendre polynomials — 
e.g., p. 175 of Ref. 2 - it follows directly that 
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0 1 2 3 4 5 

k 


FIGURE 1. 


MAPPING OF LEGENDRE INDICES k AND SL ONTO THE 
SINGLE INDEX m. tHE ORDERING OF THE MAPPING 
FOLLOWS THE INDICATED PATH. 



D/2 P/2 

r t 


(9) 


[¥ yk (y)¥ zJl (z)] 2 dydz = 1 ; 

-D/2 -P/2 

hence, the definitions of Eqs. (6) to (8) satisfy the normaliza- 
tion condition implied by Eq. (2). . For some of the work to 
follow, it will be convenient to have expressions analogous to 
Eqs. (1), (3), and (5) written directly in terms of the Legendre 
polynomials of Eqs. (6) to (8) using the double indices k,il: 


1 

FD 


P C (y ’ z) - j 0 j o a MV (y) "^ (Z) 


D/2 F/2 


kJl FD 


-D/2 -F/2 


F c (y,z)Y yk (y) 1 l' zJl (z)dydz 


D/2 F/2 


1 

FD 


00 CO 


F*(y,z)dydz 


-D/2 -F/2 


I l 

k=0 &=0 K * 


( 10 ) 

( 11 ) 

( 12 ) 


One of the main reasons for our use of Legendre polynomials 
in the series representation of FQ(y,z) is the fact that the 
two-dimensional Fourier transform of Eq. (10) can be expressed 
simply in terms of known functions. Using a caret to denote the 
Fourier transform of a function as in Ref. 1, we can express the 
Fourier transforms of the Legendre expansion functions of Eqs. 
(7) and (8) as 


¥ 


yk 


(gj ) 


A 


. 0 ° 

H' yk (y)exp(-i2iTg 1 y)dy 

— 00 


= (2k+l) 1 " 2 (F/2)P k (Fg 1 /2) 

= (-i) k F(2k+l) 1 ' s j k (TrFg 1 ) (13) 
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and 


^ S 2 ) = ^ z ^(z)exp(-i27rg 2 z)dz 


= ( 2£+l )^(D/2 ) (Dg 2 /2 ) 
= (-1 ) ^D(2£+l )^j £ (TrDg 2 ) 


where 3 for n = k and we have defined 


r oo 

P n (g) = J P n ( x ) exp ( -12irgx ) dx . 


( 14 ) 


(15) 


In arriving at Eqs. (13) and (14), we have used Eq. (1) on p. 122 
of Ref. 3 and the definition 


J n (x) = [tt/(2x) ]^J n+1/2 (x) 


( 16 ) 


where j n (x) is the spherical Bessel function of the first kind of 
order n as defined on p. 437 of Ref. 4 in terms of the Bessel 
functions Jn+l/2^ x ) of the first kind of order n+1/2. Prom 
Eqs. (1), (10;, (13)j and (14) we obtain the desired representa- 
tion of the two-dimensional Fourier transform of Fc(y,z) expressed 
in terms of the expansion coefficients defined by Eqs. (3) and 



.00 

P Q (y ,z)exp[-i2Tr(g 1 y+g 2 z) ]dydz 

_oo 



00 oo 

= FD l l (-i) K+ a [(2k+i)(2£+l)] a j (TrPgJj* (TrDgJ. 

k=0 1=0 KX/ k i £ 2 

(17) 
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Fourier-Series Coefficients of Reciprocal Mesh Stiffness and 
Load-Dependent Component of Static Transmission Error 

Equation (77) of Ref. 1 provides a general expression for 
the Pourier-series coefficients of the reciprocal mesh stiffness 
1 /K t (x) in terms of the ratios a Kn /K T as indicated. by Eq. ( 87 ) 
of Ref. 1. When tooth contact takes place over the entire 
rectangular face region 

-(P/2) < y < (P/2) 

-(D/2) < z < (D/2) 


as illustrated in the upper right-hand corner of Fig. 4 of Ref. 

1, we may use Eq. (17) to develop an expression for the Fourier- 
series coefficients of l/K T (x). In this application of the 
preceding results, we identify Fc(y,z) with the local stiffness 
of a pair of teeth per unit length of line of contact KipQ(y,z) 
expressed in the Cartesian tooth coordinates defined by 
Eqs . (15), (16), and (19) of Ref. 1. When we write the expansion 
coefficients of Eq . (11) as 

D/2 P/2 

a ia - k \ l K TC (y > z) V ly) ' 1 , z/ z)dydz ’ (18) 

-D/2 -F/2 


where Y yk ( y) and Y z& (z) are defined by Eqs. (7) and (8) respec- 
tively, we obtain from Eq. (17) the two-dimensional Fourier 
transform of K TC (y,z) defined by Eq. (82) of Ref. 1: 



(g x ,g 2 ) 


00 00 

PD l l (-i) k+£ a, [(2k+l) (2£+l)]' s x 

k=0 £=0 


j^TrFgi )0 £ UDg 2 ) . 


(19) 


Since 

j 0 (0) = 1 (20a) 

and 


j n (0) = 0, n = 1 , 2 , • • • 


(20b) 
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(Ref. 4, p. 437), we have from Eq. (19) 

K tc ( 0,0) = FDa 00 . (21) 


Therefore, from Eqs. (19), (21), and Eq. (86) of Ref. 1, It 
follows that 


00 00 

(“ K n /if T ) = l I (-i) k+t '(a k /a 00 )[(2k+l)(2Wl)] !s 
n 1 k=0 1=0 KJ6 


x J k CmrPL/(AA) ]j £ (mrL/A) . (22) 


The Fourier-series coefficients of 1/Krp(x) can be approxi- 
mated to any desired degree of accuracy by combining Eq. (22) 
with Eqs. (76), (77), and ( 87 ) of Ref. 1. For cases where M in 
Eq. (77) of Ref. 1 is taken larger than unity, we see from Eqs. 
(76), (77), and ( 87 ) that one or more discrete convolutions of 
a Kn/^T with itself are required. These convolutions generally 
will have to be carried out numerically. 

For the case where M is taken to be unity, the Fourier-series 
coefficients of 1/K rp (x) may be approximated using Eq. (77) of 
Ref. 1 by 1 


a (l/K)n “ (K T )- 1 C2(<x Kn /K Ir )t° : ' - (c^/Kj) [1] ] 


' (K T rl[2{ n,0 - ( “Kn /K T )] < 2 3) 

where Eqs. ( 87 ) and ( 76 ) of Ref. 1 have been used, and where 
5 n ,0 is Kronecker's delta. To obtain an expression for Kt for 
use in Eq. (23), we combine Eq . (21) with Eq. ( 85 ) of Ref. 1 
to obtain 


K T = a o 0 ^ EL/ A ) seci4 b 
= K TC (FL/A )sec^ b 


(24) 


where K>pc = a 00 is the average value of the local tooth-pair 
stiffness: 
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(25) 


K 


A JL 

TC ' PD 


D/2 

' 

-D/2 


P/2 

K TC (y,z)dydz 

-P/2. 
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as may be seen from Eqs. (7), (8), and (18) evaluated for k=0, 
£.=0, and from the fact that P Q (x) = 1. 

Our desired expression for the Fourier-series coefficients 
a Wn °I the l° a( ^ dependent component £y(x) 0 of the static trans- 
mission error is obtained, according to Eq. (68) of Ref. 1, by 
multiplying the total mesh loading W 0 by the Fourier-series 
coefficients of l/K T (x). According to Eqs. (22) and (23), for 
M = 1 this expression is 


a 


Wn = (W ° /K T> 26 n,0 


- I 

k=0 £,=0 


l (-i) k+£ (a 


k5, /a o o ^ 


[( 2 k+l) (25,+D ]^J k [mrFL/(AA) ]j £ (nirL/A)| (26) 

where K T is given by Eqs. (24) and (25). 

Discussion . The spherical Bessel functions in Eq. (26) 
depend on two dimensionless parameters, L/A and FL/(AA), which 
will appear repeatedly throughout the remainder of the paper. 

The parameter L/A is the profile or transverse contact ratio 
whereas FL/( AA) is the axial contact ratio [5,6], According 
to Fig. 2 of Ref. 1, L/A is the length L of the path of contact 
divided by the tooth spacing A, both measured on a line defined 
by the intersection of a plane normal to the gear axes (trans- 
verse plane [ 5 ,<?]) and the plane of contact. Furthermore, 
according to Eq. (D5) of Ref. 1, we have A = L/tani!^. If we 
denote by A a the tooth spacing in the axial direction measured 
in the plane of contact (axial pitch [ 5 ]) then from Fig. 2 of 
Ref. 1 we see that 

tamj^ = A/A a , (27) 


hence, we have 

FL/(AA) = F/ A (28) 

cl 
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which is the ratio of the face width to the axial pitch - i.e., 
the axial contact ratio (Ref. 6, p. 211). For helical gears, 
the axial contact ratio can be interpreted as the average 
number of teeth in contact across their full depth D, as illus- 
trated in the upper right-hand corner of Fig. 4 of Ref. 1. For 
spur gears, the transverse contact ratio can be interpreted as 
the average number of teeth in full contact across the gear face. 

Table 1 gives approximate percentage errors in the expres- 
sions for reciprocal mesh stiffness - Eqs. (71) to (74) of 
Ref* 1 — as a function of the truncation parameter M and the 
contact ratio. In interpreting Table 1, the contact ratio 
chosen should be the larger of L/A or FL/(AA ) . For a typical 
spur gear transverse contact ratio of (L/A) = 1.5, it may be 
seen from Table 1 that the value of M = 1 used in deriving 
Eq. (26) provides adequate accuracy for most vibration excitation 
calculations. For helical gears with larger (axial) contact 
ratios, the accuracy of Eq. (26) is better than that for the 
spur gear case. From Table 1, we see that the formulas for 
reciprocal mesh stiffness given by Eqs. (72) to (74) of Ref. 1 
all yield an error of about one percent for the contact ratios 
cited in connection with those equations. The method used to 
calculate Table 1 is described in Appendix A of this report. 

TABLE I. - APPROXIMATE PERCENTAGE ERRORS IN CALCULATION OF 
RECIPROCAL MESH STIFFNESS AS A FUNCITON OF CONTACT RATIO 

AND TRUNCATION PARAMETER M. 


Contact 

Ratio 

1 

M 

2 

3 

1.5 


11.1 

3.7 

1.23 

2.5 


4.0 

0.8 

0.16 

5.0 


1.0 

0.1 

0.01 


Application of Eq. (26) to helical gears involves a doubly 
infinite summation over the indices k and A. This summation 
arises from expansion of the local tooth-pair stiffness Kirc(y,z) 
in the two-dimensional Legendre series of Eq. (10), where the 
expansion coefficients are given by Eq. (18). From deflection 
measurements made on spur and helical gears — e.g.. Figs. 5 to 8 
of Ref. 7 and Fig. 10 of Ref. 8 (Part 2) — we would expect 
excellent accuracy in our expansion of Kijc(y>z) when the summa- 
tion over the profile (transverse) index l is truncated at £=4. 
Such a truncation will result in a least-squares fourth-degree 
polynomial approximation of the function Krp C (y,z) in the variable 
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z. For helical gears, the number of terms required in the sum- 
mation over k for a given accuracy requirement will depend on 
the ratio F/D - as this ratio increases, it can be seen from 
tooth deflection measurements (Refs. 8-10) that an increasing 
number of terms will be required. 

For application of Eq. (26) to spur gears, we note from 
Eq. (C4 ) of Ref. 1 - or from Fig. 4 of the same reference - 
that for zero helix angles we have (L/A) = 0. Thus, using 
Eqs . (20a) and (20b) we see that Eq. (26) reduces in the ease 
of spur gears to 


a 


Wn 


(w 0 /k t ) 



00 2 

1 (-i)*(a /a 00 )(2&+l)^ jl (mrL/A) 

£=0 


(29) 


which involves only a single summation over % . Using the fact 
that P (x) = 1, we see from Eqs. (7) and (18) that for spur gears 
the expansion coefficients in Eq. (29) can be expressed as 


a 


oil 



F/2 

| K TC (y,z)dy 
-F/2 


' 1 W z,dz > 


(30) 


which depends only on the "axial average" of the local tooth-pair 
stiffness K TC (y,z) — see the sketch in the upper right-hand 
corner of Fig. 4 of Ref. 1. 

We can obtain additional insight into Eq. (29) by recognizing 
that, to a first order of approximation, we can assume that 
K TC (y,z) is a constant. In this case - as may be seen from 
Eqs. (8) and (30) - Eq. (29) reduces to 

“Wn ' ( V‘V C25 n,0‘V m,L/A)] . (3D 


where (Ref. 4, p. 438), we have 


j 0 (mrL/A) 


sin(nirL/A) 

nirL/A 


(32) 


According to Eq. (32), if L/A is an integer we have j 0 (nirL/A) = 0 
for n=l,2,3,***; hence, from Eq. (31) it follows that if L/A is 
an integer, we have Wn = 0 for n = 1,2,3»*** — i.e., the tooth 
deformation contributions to the harmonics of the static trans- 
mission error at the tooth— meshing frequencies vanish. Reductions 
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in gear noise for integral contact ratios have been observed 
experimentally - e.g.. Pig. 9 of Ref. 11. However, for accurate 
spur gear calculations, terms at least through £=4 should be 
retained in the expansion of the local stiffness Ki>c(y,z); in 
addition, terms at least through M=3 should be retained in Eq. 
(77) of Ref. 1. 

In the case of helical gears when we assume that Kijn(y,z) 
is a constant, Eq. (26) reduces to 

a Wn ~ ( W o/K T ){26 nj0 -j o [nTTFL/(AA)]j 0 (n7rL/A)} , (33) 


as may be seen from Eqs . (7), (8), and (18). Prom Eqs. (32) and 
(33) j we see that if L/A is an integer, we have a w =0 for 
n = 1,2,3, in this case also; moreover, since Wn 


j 0 [mrPL/ (AA) ] 


sin[nTrFL/(AA) ] 
mrFL/( A A ) 


( 34 ) 


it follows that if FL/ (AA ) is an integer, we also have a y n = 0 
for the tooth-meshing harmonics n = 1,2,3, •••. Thus, in the 
helical gear case when Kmg(y,z) is assumed to be a constant 
we require that either L/A or FL/(AA) be an integer for the 
vanishing of the tooth deformation contribution to the tooth- 
mesing harmonics. 

These integral contact ratio criteria can be understood with 
the aid of Fig. 2, which illustrates the zone of contact, as 
shown in Fig. 2 of Ref. 1, but for cases of integral transverse 
and axial contact ratios (L/A) = 2 and (FL/AA) = 3. From Fig. 2, 
we see that if either contact ratio is an integer, there is no 
fluctuation in the total length of lines of contact within the 
zone of contact as the gears rotate. Consequently, if the local 
tooth-pair stiffness per unit length of line of contact K TC (y,z) 
is a constant, then when either the transverse or the axial 
contact ratio is a constant there is no variation in the total 
stiffness of the mesh as the gears rotate. Hence, for this case, 
the tooth deformation component of the static transmission error 
is the constant value 

a W0 = VK t . (35) 

as may be seen from Eqs. (20a) and (33). Notice that the value 
of the d-c component given by Eq. (35) also is predicted from 
the more general form of a ^ n given by Eq. (26) — as one may see 
from Eqs. (20a) and (20b). 
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FIGURE 2. BASE PLANE ZONES OF CONTACT WITH NONFLUCTUATING 
TOTAL LENGTH OF LINES OF CONTACT. ZONE OF CON- 
TACT IN UPPER SKETCH SATISFIES (L/A) = 2. ZONE 
OF CONTACT IN LOWER SKETCH SATISFIES (FL/AA) = 3 



For accurate calculations involving helical gears with axial 
contact ratios less than 5, values of M larger than unity should 
be used in Eq. (77) of Ref. 1, as indicated by Eqs. (72) to (74) 
of Ref. 1 and Table 1; moreover, additional terms beyond k=0 
and 1=0 should be used in the summations in Eq. (26), as we 
have indicated above. 

Finally, we note - e.g., from Eq. (16) and Eq. (37.4) on 
p. 139 of Ref. 12 - that as x-*-°°, the spherical Bessel functions 
j n (x) have the asymptotic form 

J m ( x ) = x ~ 1 cos[x-(m+l)ir/2] + 0(l/x 2 ), x+°°. ( 36 ) 


It follows from Eqs. (26) and ( 36 ) that for large n, the envelope 
of a Wn decays as 1/n 2 . However, for sufficiently large n this 
result will break down due to the fact that the zones of contact 
on the faces of real gears have rounded corners. 


Harmonic Analysis of Mean and Random Components of 
Static Transmission Error Caused by Deviations of 
Tooth Faces from Perfect Involute Surfaces 

Orthogonal expansions of deviations of tooth- faces 
from perfect involute surfaces 

Let us consider expansion of the stiffness-weighted mean and 
random deviations of the tooth faces from perfect involute 
surfaces : 


m 


(•) 

KC 


<y.z) = I a<° 

m m 


y 

KCm 


(y 3 z) 


and 


e 


(•) 
KC j 


(y 3 z) 


= l b 


m 


(•) 
j >m 


KCm 


(y.z) 


(37) 


(38) 


which are Eqs. (126) and (131) of Ref. 1. The mean deviations 
m KC (y> z ) and expansion coefficients are independent of 

tooth number j, whereas the random deviations e£*}(y,z) and their 

x\0 J 

expansion coefficients generally differ from tooth to tooth and 
therefore depend on j . We shall use the same expansion functions 
^KCm y * z ^ for b °th the mean and the random deviations: 

^KCm^* 2 ^ = K Tc’ 1, Cn/ y,z ^ » (39) 
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where 



A _1_ [[ 
A' 

% 


K TC (y,z)dydz 


(40) 


is the average value of the local tooth-pair stiffness and A* is 
the area of the tooth contact region . 

In the case of the mean deviations m^*^(y,z) of Eq. (37) > 

/ # \ __ AO 

if we identify m^ ; (y,z)/K TC with the function F c (y,z) in Eq. 

(1), we see from Eqs . (37) and (39) that we may apply the general 
results of Eqs. (1) to (5) to the present situation — i.e.. 





(y s z)/K TC 


'l'cm (y J z ) dydz 


(41) 


1 

A' 



r 

J - 

1 


n C 


(•) 

KC 


2 


(y,z)/K TC 


dydz = 



(42) 


Similarly, in the case of the random deviations of Eq. (38), if 
we identify 5 Z ^^TC with F cOjZ) in Eq. (1), we have 


b 


(♦) 

j 


_1_ ff 

A' 

S # 



(y,z)/K TC T Cm (y,z)dydz 


(43) 


and 


_ 1 _ 
A * 


1 1 




(y,z)/K TC 


2 

dydz 



(44) 


where in Eqs. (42) and (44) we have assumed that the functions 
^Cm^ y,z ^ Form a complete set. 
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In the discussion of Eq. (2), we pointed out that the expan- 
sion functions ^Cm^y* 2 ) must b$ dimensionless. Furthermore, 
since area A' has the dimension of length squared, it follows 
from Eq. (40) that Kijiq and KipQ(y,z) have the same dimension 
(force divided by length squared). Thus, from Eq. (48) of 
Ref. 1 and Eq. (41), we see that the expansion coefficients 
( • ) 

a m have the same dimension (length) as the nonweighted mean 
tooth deviations m£*\y,z); furthermore, from Eq. (51) of Ref. 1 
and Eq. (43), we see that the b^*^ have the same dimension 

J 5 ™ m \ 

(length) as the nonweighted random tooth deviations £ qj (y,z). 
Consequently, we see from Eq. (42) that, for each value of m, 

(a^ ^) 2 is a direct measure of the contribution of the term 

a^ ^KCm^ ,z ) (37) to the stiffness-weighted mean-square 

value of the mean tooth-face deviations m^*^(y,z). Similarly, 

/ # \ 0 

we see from Eq. (44) that (b) ' ) 2 is a direct measure of the 

j 

contribution of the term bj jJ^KCm^* 2 ^ ln Eq. ( 38 ) to the 

stiffness-weighted mean-square value of the random tooth-face 
( • ) 

deviations e :(y,z). Notice from Eqs. (48) and (51) of Ref. 1 


and Eqs. (40) to (44) above that when the local tooth-pair 
stiffness K,p^(y,z) is a constant, the values of the expansion 

coefficients a f* ^ and b^‘2, become independent of the tooth-pair 
in j ^ m 


stiffness . 


The results of Eqs. (37) to (44) - together with the material 
in Sec. II. E of Ref. 1 - provide a general set of formulas for 
the contributions to the Fourier-series coefficients of the mean 
and random components of the static transmission error caused by 
tooth-face devaitions from perfect involute surfaces. These 
results also provide formulas for the power spectrum of the 
component of the static transmission error caused by random 
tooth-face deviations from perfect involute surfaces. 

Legendre polynomial expansions of tooth- face deviations 

To apply the above results to specific sets of expansion 
functions ^Cm^y* 2 )* m = 0,1,2, •••, the functions f K m(g) and 
<P m (g) must be evaluated using Eqs. (123) and (124) of Ref. 1. 

In the case of the Legendre expansion functions defined by Eqs. 

(6) to (8), it follows from Eqs. (13) s (14), and (39) that the 
functions $K m (g) of* Eq. (124) of Ref. 1 may be expressed as 
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? Km ( S ) = (L/D)sec ^ K TC $ Cm [(L/A)g, (L/D)g] 


= K T A(-i) k+£ [(2k+l)(2jl+l)] 1 ' 2 j k (7rgFL/A) jil (7TgL) , (45) 


where Eq. (125) of Ref. 1 and Eq. (24) have been used, and where 
a mapping of the indices k and Z onto the index m, such as that 
illustrated in Pig. 1, is implied by our subscript notation. 

When Eq. (45) is combined with Eq. (123) of Ref. 1, we obtain 
a general expression for $ m (g) in terms of the Fourier-series 
coefficients ot / i/k) n. the reciprocal mesh stiffness. It was 

pointed out in Sec. IC that, when the truncation parameter M is 
taken larger than unity, the discrete convolutions required to 
compute the a (i/K)n by Eq. (22), and Eqs . (76), (77), and ( 87 ) 
of Ref. 1, generally will have to be carried out numerically. 
However, when M is taken to be unity, the Fourier-series coeffi- 
cients cc(2./K)n are gi ven by the right-hand side of Eq. (26) after 
division by W , as may be seen from Eq. (68) of Ref. 1. Combining 
the resulting 0 expression for a(i/g)n with Eq. (45), and Eq. (123) 
of Ref. 1, we obtain for $ m (g) = $k&(s)> w ^ en M = 1: 




(-i) k+£ [(2k+l)(2£+l)]^ 


j k ( tt gFL/ A ) j £ ( u gL ) - 


CO / C CO 00 , . A . 

I ] I l <-i> k + * <V*-/a 

n »=_oo\ ( k ,= 0 Z'=0 


0 0 


)[(2k’+l)(2£ , +l)] 2 x 


except 

n’=0 


j k , ( n « ttFL/AA ) j z , ( n 1 ttL/A )| j k [ ( Ag-n 1 ) ttFL/AA ] j & [ ( Ag-n 1 ) ttL/A 


(46) 


where a k , A , and a 00 are the expansion coefficients for K TC (y,z) 

given by Eq. (18), and where Eqs. (20a) and (20b) have been used. 
The number of terms required in the double summation over k’ and 
& ’ , and the approximate error involved in the truncation M=l, 
were discussed in an earlier section. 
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The function $ kA (g) plays a central role in the remainder of 
our analysis. The first term j k (7rgFL/A) j. (irgL) within the 
brackets in Eq. (46) arises from the dc component a (i/K)o 

00 

of the reciprocal mesh stiffness, whereas the terms £ (•••) 

n »=-oo 

except 
n’ = 0 

arise from the fluctuating components of the reciprocal mesh 
stiffness. See Eq. (123) of Ref. 1. 

In the case of spur gears, we have (L/A) = 0 as noted earlier. 
Thus, using Eqs. (20a) and (20b), we see that Eq. (46) reduces 
for spur gears to 


$ o£ (g) = (-i) £ (2£+l) 1 ' 5 


j £ (TTgL) 


oo / 4 oo 


l H l_ (~i) £ (a„j l »/a oo )(2£ , +l)^j^ l (n , 7TL/A) ( 


n »=_oo \ * = 0 

except 
n'=0 


I- 


j A [ (Ag-n 1 )irL/A] 


)] ■ 


/J 


* kt (g) - 0, 


k = 1,2,3 


3 i 


(47a) 

(47b) 


where a Q £ is given by Eq. (30). In view of the accuracy limita- 
tions associated with the approximation M=1 in the case of spur 
gears, Eq. (47) must be regarded as a rough approximation only. 
See Table 1. 

Fourier-series coefficients of mean deviation components 

After replacing the summation index m in Eq. (130) of Ref. 1 
by the dual indices k and i, we can express the Pourier-series 
coefficients of the mean tooth-face deviation components of the 
transmission error by 


a 


(•) 

mn 


00 


I 

k=0 


00 


l 

1=0 





(48) 
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where $j c ^(n/A) can be evaluated from Eq. (46) by setting g = n/A 
whenever the accuracy provided by taking M = 1 is adequate. 

An expression for the coefficients a^£^ in Eq. (48) can be 
obtained by combining Eqs . (6) and (41) with Eq. (48) of Ref. 1: 


(•) _ 1 


‘kJt 


FD 


D/2 F/2 

f [K TC (y,z)/K TC ] m£‘ :> (y,z)'F yk (y)'F, 7( ,(z)dydz , 


J 

-D/2 -F/2 


zl 


(49) 


where 'F . (y) and 'F „(z) are the normalized Legendre polynomials 
yK zx, / \ 

defined by Eqs. (7) and (8), and where (y,z) is the unloaded 

mean deviation of the tooth faces on gear (•) from perfect 

involute surfaces as defined by Eq. (49) of Ref. 1. 

In the case of spur gears 3 we see from Eq. (47b) that Eq. 
(48) reduces to 


a 


(•) 

mn 


X/ — u 


(50) 


where $ -(n/A) can be evaluated from Eq. (47a). 
we see that the expansion coefficients required 
be expressed as 


From Eq . ( 49 ) , 
in Eq. (50) can 


a 


(•) = 1 
oi D 



F/2 

j [K TC (y,z)/K TC ] m^‘ ^ (y,z)dy>'F z£ (z)dz 
-F/2 ) 


(51) 


where we have used the fact that P Q (x) = 1. 

Fourier -series coefficients and power spectra of 
random deviation components 


The Fourier-series coefficients of the random components of 
the transmission error can be expressed in a manner similar to 
that above by replacing the summation index m in Eq. (134) of 
Ref. 1 by the indices k and l: 


a (>) = I 
rn k=0 


I B 
£=0 


k£ 




( 52 ) 
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where $^^(n/N^*^A) can be evaluated from Eq. (46) by setting 

( # ) 

g = n/N A whenever the accuracy provided by taking M=1 is 

adequate. The expansion coefficients B^^(n) in Eq. (52) are 
given in this case by Eq. (135) of Ref. 1 after substitution 
of the indices k,£ for the index m: 




, N (,) -l , V , V 

“TTY I b) ^ exp(-i2Trnj/N l ‘ * ; ) . 
N' ' j=0 JjK 


( 53 ) 


( • ) 

An expression for the coefficients b) ^ in Eq. (53) can be 
obtained by combining Eqs . (6) and (43) with Eq. (51) of Ref. 1 

D/2 P/2 


b <-> = JL 

j ,k£ PD 


-D/2 -P/2 


[K TC (y,z)/K TC ] ^ (y j z ) Y k (y) H' z£ (z)dydz , 

(54) 


where ¥ . (y) and y 0 (z) are the normalized Legendre polynomials 

yK Z)L ( • ) 

defined by Eqs. (7) and (8), and ei. (y,z) is the random com- 

o j 

ponent of the deviations of the tooth faces of gear (•) from 
perfect involute surfaces as defined by Eq. (50) of Ref. 1. 

In the ease of spur gears 3 Eqs. (52) and (54) reduce to 


a 


= l B^(n) $ nP (n/N (,) A) 
rn o J6 T o&' 


(55) 


and 


b (,) 

b j 


D/2 ( 

P/2 

1 

n 

f 

D 

)f 

J 

-D/2 ^ 

-P/2 

from 

Eqs . 

(47b) 


(56) 


which follow from Eqs. (47b) and (7), and the fact that P (x) = 1 
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Formulas for the power spectra of the random components of 
the transmission error for a given gear are given by Eqs. (137) 
to (140) of Ref. 1. Application of the formulas to the present 
situation implies a mapping of k and i onto m such as that 
illustrated in Fig. 1. These formulas can be evaluated from 

the set of coefficients b(*^. 0 given by Eq. (5*0 » or Eq. (56) for 

spur gears, and the functions <J>, 0 (g) given by Eq. (46) and 

( . ) 

evaluated at g = n/N A. 


Discussion 

For a generic pair of meshing gears there are three compon- 
ents to the deterministic contribution of the static transmission 
error — the tooth deformation component and the component from 
each of the two gears caused by the mean deviation of its tooth 
faces from perfect involute surfaces as indicated by Eq. (6l) 
of Ref. 1. Each of these three components is a periodic function 
of x with period A. Consequently, each Fourier-series coeffi- 
cient of the total deterministic contribution to the static 
transmission error is given by the sum of the Fourier-series 
coefficients from each of the three above-mentioned deterministic 
components as indicated by Eq. (65) of Ref. 1: 


a 


mn 


= a 


Wn 


+ a 


( 1 ) 

'mn 


+ a 


( 2 ) 


mn 


(57) 


Expressions for a and a.^\ (•) = (1) or (2), are given, e.g., 

wn mn 

by Eqs. (26) and (48). 

In addition, for each pair of meshing gears there are two 
components to the random contribution of the transmission error 
as indicated by Eq. (62) of Ref. 1. The Fourier-series coeffi- 
( • ) 

cients of each of these random components are given by 

Eq. (52) . The harmonic components of each of these two random 
components generally occur at different frequencies — e.g., 
when the numbers of teeth on the two meshing gears are prime 
to one another, the random components of the transmission error 
from each of two meshing gears have no nonzero harmonics that 
occur at the same frequencies [I]. 

The rates of convergence of the series in Eqs. (48) and (52) 
are governed, respectively, by the rates of convergence of the 
expansions in Eqs. (37) and (38) when the expansion functions 
^Cm^y* 2 ) in Eq. (39) are chosen to be the Legendre polynomials 
of Eqs. (6)_ to (8). Dividing Eqs. (37) and (38) by the mean 
stiffness Kipc> then using Eqs. (6) to (8), (39), and Eqs. (48) 
and (51) of Ref. 1, we obtain 
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( 58 ) 



( y j z ) /K|j,q 



(y>z) = 


l l a^* ) [(2k+l)(2Jl+l)] 1 ' 5 P (2y/F)P.(2z/D) 
k=0 £=0 KJC k x, 

[K TC (y,z)/K TC ]e^j^ (y,z) = 

00 °o 

1 l b i k o[(2k+l)(2£+l)] J5 P (2y/F)P (2z/D) , (59) 

k=0 £=0 J 5 • ~ 


where we have used double subscripts k£ Instead of single sub- 
scripts m, as before, and where the expansion coefficients are 
given by Eqs. (49) and (54). The first three Legendre polynomials 
are [2] 

P 0 (x) = 1 


Pj(x) = x 

P 2 (x) = |(3x 2 -1). (60) 


For the purpose of interpreting the various low-order terms 
in Eqs. (58) and (59), if we assume that the local tooth-pair 
stiffness is a constant — i.e., KrpQ(y,z) = Kijq — then we may 
identify these low-order terms with qualitative "error" descrip- 
tions used in the gear-engineering literature [6,13]. In this 
interpretation, the left-hand sides of Eqs. (58) and (59) repre- 
sent, respectively, the mean and random deviations of the tooth 
faces from perfect involute surfaces. When we assume that 
K TC (y,z) = Ktc > the term k = 0, i = 0 in Eq. (58) represents a 
constant displacement of every tooth on a given gear — i.e., a 
rigid body rotation of the gear. However, since the terms in 
Eq. (59) depend on tooth number j, the term k = 0, l = 0 in Eq. 
(59) represents the pitch error of the jth tooth — i.e., the 
tooth-spacing error. According to the sketch in the upper right- 
hand corner of Fig. 4 of Ref. 1, the index k in Eqs. (58) and 
(59) is the lead "error" index; whereas the index Z is the 
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profile "error" index. Consequently, the terms k = 0, & = 1 in 
Eqs. (58) and (59) represent involute slope deviations from 
perfect involute surfaces, since in this case the right-hand 
sides of Eqs. ( 58 ) and (59) are linear functions of the variable 
z illustrated in Fig. 4 of Ref. 1. Direct extension of this line 
of reasoning yields the following interpretations, which are 
valid when we assume that K^q(i/ 3 z) = K^q: 


k = 0 3 Z = 0 
k = 0 3 1 = 1 
k = l 3 1 = 0 
k = o, z = 2 
k = 1 3 1 = 1 

k = 2 3 Z = 0 


tooth-spacing deviations 

pure involute slope deviations 

pure lead mismatch deviations 

pure involute hollow (fullness ) deviations 

combined lead mismatch-involute slope 
deviations 

pure lead hollow (crowning) deviations . 


From the above interpretations and their extensions to the higher- 
order Legendre polynomials we may estimate the number of terms 
required for adequate engineering approximations in Eqs. (58) 
and (59), and hence the number of terms required in Eqs. (48) 
and (52) and in their respective spur gear counterparts, Eqs. 

(50) and (55). Furthermore, from the truncated sums of the 
squares of the expansion coefficients, we can use Eqs. (42) and 
(44) to determine the mean-square errors involved in truncations 
of the expansions in Eqs. ( 58 ) and (59). 

Each term in the summations of Eqs. (48) and (52) contains 
two factors — an expansion coefficient and the function 
<j>k 5 ,(§)* The expansion coefficients characterize the stiffness- 
weighted derivations of the tooth faces from perfect involute 
surfaces, whereas the function $k£(g) is independent of these 
deviations. Since these stiffness-weighted deviations [together 
with the load-dependent component of Eq. (26)] determine the 
vibratory excitation, we may regard each term in Eqs. (48) and 
(52) as a kind of transfer relationship that governs the vibra- 
tory excitation — the expansion coefficients characterize the 
"inputs" and the function $k£,(g), which is common to both Eq. 

(48) and Eq. (52), characterizes the effect of the meshing 
action of the gear pair on the inputs. ^The product of one of 
these expansion coefficient inputs and <f>k£,(g) yields the 



mesh-attenuated input. Similar conclusions can be drawn from the 
power-spectrum representation of the random component described 
by Eq. (140) of Ref. 1. 

It is appropriate to use the term ^"mesh transfer function" 
to describe each of the various terms <f>k£,(s). These mesh transfer 
functions depend only on the local tooth-pair stiffness and the 
nominal design parameters of the gear mesh, as may be seen from 
Eq. (46) for the case where M=1 . This transfer function interpre- 
tation of $ks,(g) is further supported by the fact that every 
$k$,(g) is dimensionless. Further discussion of the mesh transfer 
functions <J>k 5 ,(g) is provided in the next major section of this 
report . 

Alternative expansion method of tooth-face deviations 
from perfect involute surfaces 

We shall see later that from a certain theoretical point of 
view with important practical implications a set of expansion 
functions ^KCm^y* 2 ) different from that indicated by Eq. (39) is 
preferable for representing the deviations of tooth faces from 
perfect involute surfaces. Let us consider now, instead of the 
set of expansion functions in Eq. (39), the set defined by 

^KCm^* 2 ^ = K TC^ y,z ^ ^Cm^ y J z ^ » (6l) 


for use in representing the stiffness-weighted mean and random 
tooth-face deviations by Eqs. (37) and (38) respectively. The 
set of functions 4 , Q m (y,z) Eqs. (39) and ( 6 l) are assumed to 
be the same; however, in Eq. (39) we have weighted the expansion 
functions with the average tooth-pair stiffness defined by Eq. 

(40), whereas in Eq. ( 6 l) this average stiffness has been replaced 
by the local tooth-pair stiffness which is a function of the tooth- 
face coordinates y,z. Notice that when K Tr ,(y,z) is a constant 
value, Eq. ( 6 l) becomes identical with Eq. (39). 

By comparing Eqs. (37) and ( 38 ) with Eqs. (48) and (51) of 
Ref. 1 , we see that when the expansion functions 4 , KQ m (y,z) are 
defined as in Eq. ( 6 l) the expansions of Eqs. (37) and (38) are 
equivalent to expansions of the nonstiffness-weighted mean and 
random deviations directly in terms of the functions (y,z) - 
i.e., 



(y,z) = l a 
m 


(•) 

m 


' f Cm (y > z) 


( 62 ) 
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and 


e 0j )(y ’ z) = J b j ]m ' f 0m (y * z) • (63) 

o m 0 3 

We may again apply the general results of Eqs. (1) to (5) to the 

expansions of Eqs. (62) and (63). By identifying F P (y,z) with 
/ \ / \ ^ 

m^ ; (y,z) and e^. (y,z) in succession, we see from Eqs. (3) and 
(5) that our expansion coefficients for the mean deviation 
mQ * ^ ( y , z ) are 


a 


(•) 

m 




(y 5 z) 'J' Cm (y,z)dydz , 


(64) 


which satisfy 


1 

A 1 


[ [m£ * ) (y,z) ] 2 dydz = £ (a£*^) , 

• m 


(65) 


^ , 


whereas our expansion coefficients for the random deviations 
£q j ^ (y > z ) are 


.(O = 

J,m ** 


1_ 

A T 


45 } (y,z) ^ Cm (y j z ) dydz , 


( 66 ) 




which satisfy 


1 

A’ 


[ £ ci ^ (y > z ) 3 2 dyd 


z = 






(67) 


When the expansion functions 't'cm^> z ^ are chosen to be the 
Legendre polynomials as in Eqs. (6) and (8), the expansion 
coefficients of Eqs. (64) and (66) take the form of Eq. (11) 

when F c (y,z) is identified with m^*^(y,z) and e£j^(y,z), 

respectively . 
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In addition to the expansion coefficients a^* and bj’]., 
the functions $ m (g) are required for evaluation of the 
Fourier-series coefficients of the mean and random deviation 
components of the tooth faces, and for the power spectrum of 
the random components as may be seen from Eqs. (130), (13*1), 
and (1*10) of Ref. 1 . From Eqs. (123), (12*1), and (125) of 
Ref. 1 , we see that to evaluate $ m (g), we require the two- 
dimensional Fourier transform Y KCm (g. ,g.) of V KCm (y,z), which 
in the case of Eq. ( 6 l) is expressed^ 


^KCm (s i ,s 2 ) 


J j K TC (y ’ z) ' 1 , Cm^ y,z ^ 


x exp[-i 2 ir(g 1 y+g 2 z) ]dydz . ( 68 ) 


Thus, when ^KCm^y* 2 ) chosen to take the form of Eq. ( 6 l), 
^KCm^S ! j S 2 ) depends on the functional behavior of the local 
tooth-pair stiffness K<r C (y,z) and therefore cannot be expressed 
in closed form for a generic tooth-pair stiffness function. In 
contrast, when ^KCm^* 2 ) was chosen to take the form of Eq. (39), 
^KCm^Si>g 2 ) could be expressed using Eqs. (13) and ( 1 * 1 ) as a 
function proportional to the product of two spherical Bessel 
functions [when TQ m (y,z) was chosen to have the form of Eqs. 

( 6 ) to ( 8 )]. This simplification resulted in the relatively 
simple expression for f Km (g) given by Eq. (45) that, in turn, 
permitted us to obtain the general expression for $v^(g) given 
by Eq. (46). However, for specific tooth-pair stiffness 
functions KmQ(y,z), use of the expansion functions Tj{Q m (y,z) 
defined by Eq.^( 6 l) is practical when it is anticipated that 
the functions $ m (g) are to be evaluated numerically at the 

required values of g = n/N^A, n = 0 , ± 1 , ± 2 ,***. 


Constant Tooth-Pair Stiffness Formulas 

In this section, we summarize the special cases of the 
formulas derived in the preceding sections that apply to situa- 
tions where the local tooth-pair stiffness per unit length of 
line of contact is assumed to be a constant — i.e., Kij 0 (y,z) = 
KipQ, where is defined by Eq. (40). For these situations, 
the stiffness-weighted expansion functions given by Eqs. (39) 
and ( 6 l) are identical, and the alternative expansion method 
outlined in the preceding section reduces identically to the 
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method described earlier in the paper. In all of the considera- 
tions in this section, it is assumed that tooth contact takes 
place over the rectangular tooth-face region shown in the upper 
right-hand corner of Fig. 4 of Ref. 1. Formulas for spur gears 
are reduced in a separate section from those derived for helical 
gears . 

Constant tooth-pair stiffness formulas for helical gears 

Fourier-series coefficients of mesh stiffness and load- 
dependent component . Equations (68), (77), and ( 87 ) of Ref. 1 
provide a set of formulas for the Fourier-series coefficients of 
the load-dependent component ?y(x) 0 of the static transmission 
error in terms of the Fourier-series coefficients aj( n of the 
instantaneous total_mesh stiffness Kip(x) divided by the time- 
averaged stiffness Kt* For the case where Kijig(y,z) = we 

see from Eqs . (22), (32), and (3*0 that these latter normalized 
Fourier-series coefficients reduce to 


, /y \ _ sin(nirFL/AA) sin(nirL/A) 
K Kn 7 T 7 nuFL/AA nuL/A 


(69) 


where the vanishing of all coefficients ay^ in Eq. (22) except 
a Q0 follows directly from the orthogonality of the two-dimensional 
Legendre expansion functions H'yk(y ) ^ z s, ( z ) given by Eqs. (6) to 
(8), the fact that 'kyQ (yJ'i'^o = 1> and the expression for a^£ 
given by Eq. ( 18 ) with K TC ty,z) = K^q . For the case where 
M=1 in our approximation of the reciprocal mesh stiffness given 
by Eq. (71) of Ref. 1, it follows from Eqs. (23) and (69), and 
Eq. (68) of Ref. 1 that when KrpQ(y,z) = Kij>q , the Fourier-series 
coefficients of the load-dependent component £ w (x) 0 can be 
expressed as 


a Wn ~ ( W 


26 


n,0 


sin(mrFL/AA) sin(nirL/A) 


nirFL/AA 


niTL/A 


(70) 


hence, since (sinx/x) = 0 at x = 0, Eq. (70) is equivalent to 


a Wo ~ W o //K T 


(71a) 


a 


Wn 


- < W ^nlT/A 7 ^ * » 


(71b) 
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In the above formulas, W 0 is the total (constant) force trans- 
mitted by the mesh as illustrated in Fig. 1 of Ref. 1, and Kij is 
the time-averaged value of the total mesh stiffness defined 
by Eq. (56) of Ref. 1, which is related to the_ tooth-pair 
stiffness per unit length of line of contact Kip C by Eq. (24). 

Mesh transfer functions. When K-rc(y,z) = K TC , all expansion 
coefficients a k £ of K TC (y,z) given by Eq. (18) vanish except 
a 0 o as^indicated above. For this case, the mesh transfer func- 
tions <f> k £(g) given by Eq. (46) for M=1 reduce to 

$ k£ ( S } = (-i) k+£ [(2k+l)(2il+l)] i5 {j k (7rgFL/A)j £ (7rgL) 


00 

I j 0 ( n ' ttFL/AA ) j (n ’ ttL/A ) J [ ( Ag-n ' ) ttFL/AA ] j [ ( Ag-n ’ ) ttL/A ] } . 

n f = —00 ^ ^ 

except 

n*=0 ( 72 ) 


Fourier-s eries coefficients of mean deviation components . 
Equation (48) expresses the Fourier-series coefficients of the 

mean deviation transmission error component from gear (•) 

in terms of the expansion coefficients a k £^ and the mesh transfer 
functions $ k £(n/A). In situations where KipQ(y,z) = K<pc, we see 

from Eq. ( 58 ) that the a k £^ are the expansion coefficients of 

the mean deviation mQ*^(y,z) of the tooth faces of gear (•) from 
perfect involute surfaces (without any stiffness weighting) as 
defined by Eq. (49) of Ref. 1 — i.e.. 


. n w 00 . 

mi* ; (y,z) = l l a^* ) [(2k+l)(2iL+l)] is P (2y/F)P (2z/D) , (73) 

0 k=0 1=0 K £ 


where, from Eqs . (7), ( 8 ), and (49), we see that the expansion 
coefficients are obtained from mQ*^(y,z) by 


.(•) = [ ( 2 k+l ) ( 2 &+ 1 ) 
kH FD 



F/2 

j m ( l* ) (y,z)P k (2y/F)P £ (2z/D)dydz 
-F/2 

(74) 
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where Pj c (x) and P&(x) are the Legendre polynomials as defined in 
Ref. 2 and Eq. (60). Prom Eq. (48) of Ref. 1 and Eq. (42), we 

see that when Kmp(y,z) = Km P , the expansion coefficients a/ * ^ of 
Eq. (74) satisfy kZ 


D/2 P/2 

a 1 l 

-D/2 -P/2 


[m£ * ^ (y,z) ] 2 dydz 


°° °° ( • ^ 2 
■ l l (a M> 


k=0 a=Q 


(75) 


Therefore, the squares of the expansion coefficients a£j^ are 
direct measures of the contributions of the various terms in 

Eq. (73) to the mean-square value of mp*^(y,z) given by the 
left-hand side of Eq. (75). 0 

In addition to the expansion coefficients given by Eq. (74), 
our expression, Eq. (48), for the Pourier-series coefficients of 
the mean deviation component requires the mesh transfer functions 
$k£,(s) evaluated at g = n/A — which we find from Eq . (72) to be 
for the case M=1 : 

«|> kJl (n/A) = (-i) K X '[(2k+1)(2£+1)] 2 {j k (n7rFL/AA)j £ (mrL/A) 


l J 0 ( n * ttFL/A A ) j 0 ( n • ttL/A ) j . [ ( n-n • ) ttPL/A A ] j . [ ( n-n » ) ttL/A ] } . 
n 1 =-“ 
except 
n * =0 


(76) 


Let us consider <j> oo (n/A): 


£ _ sin (nuFL/ AA ) sin(nirL/A) 

J nuFL/AA nirL/A 


0 0 
00 

- I 

n * =-°° 

except 

n’=0 


n'lrPL/AA n'irL/A (n-n 1 )ttFL/AA (n-n*) ttL/A 


sin (n * ttFL/ AA ) sin(n 'ttL/A ) sin[ (n-n 1 ) ttFL/A A] sin[ (n-n 1 ) ttL/A] , 

A 

(77) 
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where we have used the fact that j 0 (x) = (sinx)/x as noted 
earlier. In Appendix B, it is shown that the most important 
term within the summation in Eq. (77) is the term n'=n. 
Retaining only this one "correction term" gives 


<f>oc>/A) = 


1 , 


0 , 


n = 0 


n ? 0 


( 78 ) 


since at x = 0, we have (sinx)/x = 1. 

Although the derivation of Eqs. (77) and (78) involved sev- 
eral approximations, the result given by Eq. (78) is exact. To 
verify the validity of Eq. (78), consider the harmonic n = 0 
first. Since P 0 (x) = 1, we see from Eq. (73) that the contri- 
bution to rriQ ^ (y,z) arising from the term aoo^ is simply a[ o ^ — 
i.e., the same constant value added uniformly to the face of 
every tooth on gear (•). To interpret this result, we have 
from Eqs. (9) to (11) of Ref. 1 or equivalently from summing 
Eq . (8c) of Ref. 1 over all tooth pairs j in contact: 


W = l W (x) 

j J 



(x,y)[c(x)- n 5 1) (x,y)-n? 2) (x,y)]dy . 
J J 


(79) 


For either of the gears (1) or (2) in contact, the tooth surface 
( • ) 

deviations n. are the superposition of mean and random compon- 

J 


ents as may be seen from Eq. 

( • \ 

tion of a mean component ao 0 


(50) of Ref. 1. Hence, 

( • ) 

to the deviations n : 

J 


the addi- 
of every 


tooth j on either gear (1) or gear (2) must give rise to the 

addition of the same constant value aoo^ to the transmission 
error ?(x) within the ‘brackets in Eq. (79) if the load W is to 

remain unchanged. But the addition of a constant value a^o^ to 

?(x) results in an increment of the same value aoo^ to the dc 
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( • ) 

Fourier-series coefficient a^ 0 as may be seen from Eq. (6 7) of 
Ref. 1 when n = 0. Examination of this same dc contribution in 

Eq. (48) from the term aoo^ gives aoo^ = aoo^$ 0 oC°) after sub- 
stitution of ao o ^ for in the left-hand side. We therefore 

must have $ 00 (O) = 1 as indicated by Eq. (78). Furthermore, we 


have seen from Eq 
( • ) 

ao o to the deviations 


(79) that the addition of a constant value 
( • ) 

j of every tooth on a gear is equi- 
valent to adding the same constant value to the transmission 
error c(x); hence, the only nonzero Fourier-series coefficient 


( • ) 

arising from a 00 


( . ) ( . ) 

is the dc term a„/ = aoo 


- - n * mo — a. o o • Thus, from Eq. 
(48), we must have <j> (n/A) = 0 for n ? 0 - also as indicated 
by Eq. (78). 


Conversely, we see that adding a constant value to C(x) 
within the brackets in Eq. (79) is equivalent to adding the same 

constant value to either nj 1 ^(x,y) on nj 2 ' ) (x,y) for every value 

of j if the load W is to remain unchanged. It therefore follows 
from the orthogonality of our expansion functions of the mean 

deviation components mQ*^(y,z) in Eq. (73) that the only tooth 
deviation components that can give rise to the dc component 

a^ o ^ of the transmission error are the constant deviations a|^ 


of the teeth on gears (1) and (2). • Hence, we must conclude from 
Eq. (48) that 


; M (o) = o 


whenever k and Z are not both zero. 


( 80 ) 


We shall see below that the above conclusions lead to import- 
ant consequences with regard to "random" tooth-spacing errors. 

The results indicated by Eqs . (78) and (80) are valid for arbitr- 
ary tooth-pair stiffnesses per unit length of line of contact. 
These results — which separate pure tooth-spacing errors from 
other classes of errors — constitute the main reason why the 
expansion functions of Eq. (61) are theoretically preferable to 
those of Eq. (39) in cases where K TC (y,z) is not equal to the 
constant value Kijq . The expansion functions of Eq. (6l) are 
preferable in these cases because these expansion functions 
lead to the expansions, Eqs. (62) and (63), of the tooth-face 
deviations where the first terms m = 0 always can be made to 
represent pure tooth-spacing "errors". The expansion functions 
of Eq. (39) do not possess this property for arbitrary tooth- 
stiffness functions Kijc:(y>z)j and as a result, for tooth-stif fness 
functions other than the constant value Krpc(y,z) = Kijq, the 
properties of Eqs. (78) and (80) generally will not hold when 
the expansion functions of Eq. (39) are used. 
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Fourier-series coefficients and power spectra of random 
deviation components . Equation (52) expresses the Fourier-series 
coefficients of the random deviation transmission error component 

Cr ^ from S ear (*) in terms of the mesh transfer functions 

and the discrete Fourier transforms B^^(n) of the 
expansion coefficients as defined by Eq. (53). In situa- 

tions where K TC (y,z) = K TC , we see from Eq. (59) that the 
b j ki, are the ex P anslon coefficients of the random deviations 

e Qj^(yj z ) of the faces of gear (•) from perfect involute surfaces 
as defined by Eq. (50) of Ref. 1 - i.e.. 


(y»z) = l l bj *' £ [(2k+l) (2£+l)] Js P k (2y/F)P £ (2z/D) . 
k — 0 9 j — 0 3 

( 81 ) 

From Eqs . (7), (8), and (5*0, we may express the expansion co- 
coefficients in Eq. (8l) in terms of e^(y,z) by 

o J 


b 


(•) 

j,kS, 


[(2k+l)(2 SL+l)f 2 
FD 


D/2 F/2 

1 1 

-D/2 -F/2 


e^j ) (y,z)P k (2y/F)P Jl (2z/D)dydz . 

( 82 ) 


Furthermore, from Eq. (51) of Ref. 1 and Eq. (4*1), we see that 
when K TC (y,z) = K TC , the expansion coefficients bj ^ of Eq. (82) 
satisfy 


1 

FD 


D/2 F/2 


-D/2 -F/2 


[Oq j ^ (y , z ) ] 2 dyds 


= l 


k=0 1=0 


? 0;m )2 


(83) 


Equations (137) to (140) of Ref. 1 express the discrete 
power spectrum |«rn^l 2 °f the random deviation transmission 
error component from gear (•) in terms of the expansion coeffi- 
cients and the mesh transfer functions $ kjl (n/N^ * ^A) - 
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where here a mapping of the index m onto the indices k and £ such 
as that illustrated by Pig. 1 is implied. For a specific error 
component designated by the index m or the pair of indices k£, 
the tooth-face random errors of gear ( • ) are characterized for 

( • ) 2 

this application by their mean-square values (b. ' ) and auto- 

C . ) k£ 

correlation coefficients ( j ) or discrete autospectra 
s£^(n) as indicated by Eqs. (137) to (139) of Ref. 1: 


< b M )>2 4 


, N (-, -l , , 2 

-jTy S (b j 

N' ‘ j =0 J 


(84) 


Om’o) 4 


( . ) 2 

< b M > 


- 1 


N 




N ( ‘^-l 

I 

j=0 


b ( * } b ( ‘> 

j ’ ,k£ j ’+j ,k£ 5 


(85) 


and 


, , . , • * ^-1 , , 

S M (n) = j £ 0 p ki (J ) exp 


/-i 27 rnj\ 

V7^7 ’ 


where N 


(*) 


( 86 ) 


designates the number of teeth on gear (•). The con- 
tribution |a v ; I „ of the single error component designated by the 
pair of indices k£ to the power spectrum |a^ >n '| is then given by 


ia (,) |* 4 (b <*V s { '\n) 

rn 'k£ k£ k£ 


$ kjl (n/N ( * } A) 


(87) 


which arises from a single term m’ = m" in the double summation 
in Eq. (140) of Ref. 1. When several random error components 
k£ are present, cross-spectrum .terms must be included to obtain 

all contributions to I a ^ ^ I 2 as indicated by the double summation 

in Eq. (140) of Ref. 1 — unless these error components are un- 
correlated in which case the cross-spectrum terms are zero. 
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Notice that Eq. ( 87 ) has the form of a power spectrum input- 
output relationship - e.g., Eq. (2.3*0 on p. 71 of Ref. 14. 

The input characterization is the mean-square error (b/^) 2 

k£ 

times the normalized tooth-error spectrum S^*^(n) - see Eqs. (84) 

to (87) - and the "system function" is the squared magnitude 

|4> k $/ n//N ^ ^ A l 2 of the mesh transfer function, which is independent 

of the tooth errors. The more general expression for |a^*^| 2 

given by Eq. (140) of Ref. 1 has the form of the power spectrum 
of the output of a linear system to multiple random inputs - 
e.g., Eq. (5.42) on p. 149 of Ref. 15. 

Modified forms of mesh transfer ^functions . Let us now 
consider the mesh transfer function $ k£ (g) evaluated from Eq. (72) 

at values of g = n/N^ ^A: 


$ M (n/ N C>A> . 

- n ’) ! ^] J t[(7Ty -n'j 


- I j 0 (n ’ ttFL/ AA )j (n ’ttL/ A ) j 


n 1 =-°° 
except 
n ’ =0 


ttL 

A 


( 88 ) 


At the tooth-meshing harmonics, which occur at integral values of 

P = n/N^ \ it was pointed out earlier that the exact value of 
$ 00 (P / A) is zero except for p = 0 where $oo(0) = 1, and further- 
more that $oo can be modified to yield this behavior by retaining 
only one "correction term" in the summation over n*. For the 
rotational harmonics represented by Eq. (88), this modification 
can take the form 


o ( n/N ^ * ^A) s j 


/ n ttFL \ . / n ttL \ 
AA j J o(^T A j 


- I |rect [ ( n/N^ * ^ )-n 1 ] j (n * ttFL/AA )j . (n 1 ttL/A ) 
n»=-o o ( 0 0 

except 
n ’ =0 


X J 


/ n . ,\ tt FL 1 |7 n ,\ ttL 


(89) 
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where rect x is the rectangular function [p. 29 of Ref. 16] 
defined as 


rect x = 


Prom Eq. (90) we see that for any rotational harmonic n, only 

one term in the summation in Eq. ( 89 ) can contribute — the term 
. ( . ) . 

for which | (n/N )-n'| < 1/2. Moreover, when n falls in the 

range |n/N <_ 1/2, none of the terms in the summation in 

Eq. ( 89 ) contribute. It follows, therefore, that when Eq. ( 89 ) 
is evaluated at the tooth-meshing harmonics, which occur at 

( • ) 

integer values of n/N , it satisfies exactly the condition of 
Eq. (78). 

The right-hand side of Eq. ( 89 ) will display (small) 
discontinuities at values of n in the neighborhoods of 

[(n/N^ ))-n'] = 1/2 for all n*. These discontinuities can be 

removed by replacing rect [ (n/N^ * ^ )-n ' ] in Eq. ( 89 ) by a function 
that smoothly increases and decreases the weighting of each 
term. A function ideally suited for this purpose is 



w 00 [(n/N ( ‘ ) )-n'] = rect {|-[ (n/N^ * ^ )-n 1 ] } (1+cos { ir [ (n/N^ * ^ )-n 1 ] } ) . 

(91) 


When rect[(n/N V-n’] is replaced in Eq. ( 89 ) by Eq. (91), we 

( • ) 

see that whenever n/N is not an integer, we generally will 
have two terms in the summation over n T providing nonzero con- 
tributions. However, the condition of Eq. ( 78 ) still will be 
satisfied exactly. 

Let us turn now to the general case where k and l are not 
both zero. Here, except in the vicinity of g = 0 where Eq. (80) 
must be satisfied, we have no reason for not retaining all of 
the terms in Eq. (72). The requirement of Eq. (80) can be 
satisfied while avoiding discontinuous behavior by multiplying 
the summation in Eq. (72) by 
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|Ag| < 1 


|Ag| > 1 , (92) 

which increases smoothly from zero at g = 0 to unity at |g| = 1/A, 
which is the position of the first tooth-meshing harmonic where 
no modification to the summation in Eq. (72) is sought. For 
values of |g| >_ 1/A, Eq. (92) provides no modification to the 
summation in Eq. (72). 

We may now collect the above modifications of Eq. (72) into 
a general formula that is applicable to both the mean and the 
random deviation components and that satisfies the requirements 

of Eqs . (78) and (80). Setting n/N^*^ equal to the product Ag 
in Eq. (91) gives 

w Q Q ( Ag-n ' ) = rect[|- (Ag-n’) ]{l+cos[Tr(Ag-n' )]} . (93) 


rect (Ag/2) [1-cos ( uAg) ] , 


If we now define the general weighting function 


! w o o (Ag-n * ) , 
w( Ag) , 


k = 0, l = 0 


otherwise , 


( 94 ) 


where w(Ag) and w 00 (Ag-n’) are defined by Eqs. (92) and (93) 
respectively, we then have for our general modification of 
Eq. (72): 


/s. , b-4- 0 1 A 

<f> k Z (S) ~ (-D [(2k+l)(2A+l)] 2 (j k (iTgFL/A)j Jl (TrgL) 

00 

1 {w. ? (Ag,n')j 0 (n»TTFL/AA)j 0 (n»TTL/A) 
n’=-~ KJ6 

except 
n » =0 


x J k ^ (Ag-n* )ttFL/AA]j £ [ (Ag-n ' ) ttL/A ] } I , 


( 95 ) 
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which may be used to evaluate <j> k £ in Eqs . (48), (52), and (87). 

Constant tooth-pair stiffness formulas for spur gears 

Since helical gears become spur gears when the helix angle 
is zero, the above formulas can be reduced to the case of spur 
gears by setting (L/A) = 0 as indicated by Eq. (C4) of Ref. 1. 
The resulting spur gear formulas — which are applicable to the 
case Kipc(y,z) = Kij«q ^ — are summarized below. Formulas requiring 
transfer functions $ k ^(g) must be regarded as rough approxima- 
tions only when the approximations given below for the $ k ^(g) 
are used, since these formulas have been evaluated from the 
approximation M = 1 in Eq. (71) of Ref. 1 — see Table 1. 

Fourier-series ooeffioients of mesh stiffness and load- 
dependent component. For use in Eqs. (68), (77) 5 and (87) of 
Ref. 1 for calculating the Fourier-series coefficients of the 
load-dependent component of the static transmission error, we 
have for spur gears instead of Eq. (69): 


^ a Kn /K T^ 


sin(nirL/A) 

nirL/A 


(96) 


For the case where M = 1 in our approximation of the reciprocal 
mesh stiffness given by Eq. (71) of Ref. 1, we may obtain from 
Eqs. (70) and (71) a rough approximation for the Fourier-series 
coefficients of the load-dependent component of the static trans- 
mission - error applicable to the case of spur gears: 


a Wn = <W 


26 


sin(n7TL/A ) 


n,0 nirL/A 


(97) 


or 


a Wo ~ ^o^ K t 


a Wn - -(* 0 /V Sln n ( ^7/ - Al • 


(98a) 

(98b) 


Mesh transfer functions . For spur gears, the mesh transfer 
function of Eq. (72) reduces to 
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$ o£ (g) = (-ir(2A+l)'Vj „(7rgL) - 


I J 0 (n , 7TL/A)j^[(Ag-n' )TrL/A], 


n ’=-« 
except 
n ’ = 0 


(99a) 


and 


$ k£ (g) =0^/0. 


(99b) 


In view of the fact that the above result was obtained using the 
approximation M = 1, the above formula should be used for rough 
calculations only. Furthermore, the reduction of Eq. (95) to 
the case of spur gears [given below by Eq. (110)] will give 
results generally superior to those provided by Eq. (99). 

Fourier-series coefficients of mean deviation components . 
Equation (50) expresses the Fourier-series coefficients of the 
mean deviation component of the transmission error from gear (•) 

in terms of the expansion coefficients a^ and the mesh transfer 
function cf> £ (n/A). From Fig. 4 of Ref. 1 and Eqs. (73) and (74), 
we see that these expansion coefficients actually characterize 

the axial average mQ*^(z) of the mean deviation mQ*^(y,z) of the 
tooth faces on gear (•) from perfect involute surfaces with zero 
helix angles — i.e., 

F/2 

=| j m£*^(y,z)dy . (100) 

-F/2 


To show this, we divide Eq. (73) by F, integrate over y between 
-F/2 and F/2, and then use the orthogonality of the Legendre 
polynomials and the fact that P 0 (2y/F) = 1 - thereby yielding 


* 



W - 

l a 0 ( * ) (2£+l) is P ? (2z/D) . 
£=0 A 


( 101 ) 


From Eqs. (74) and (100) and the fact that P (2y/F) = 1, we then 
have 0 
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(•) _ ( 2 &+ 1 ) h 
a oi ~ D 


I 


-D/2 


m^*^(z) P^(2z/D)dz . 


( 102 ) 


Finally, by squaring both sides of Eq. (101), then dividing the 
resulting expression by D and integrating over z between -D/2 
and D/2, and using the fact that 


( 21 + 1 ) 

D 


D/2 


-D/2 


[P p (2z/D) ] 2 dz = 1 , 


(103) 


(•) 


we see that the expansion coefficients a^ must satisfy 


1 

D 


D/2 


-D/2 


[m£*)(z)] 2 dz 


i kV ) 2 

Z=0 


(104) 


Equations (50) and (100) to (104) show that the axial average 

m£*)(z) of the mean deviation m^,*^(y,z) of the tooth faces on 

gear (•) is the appropriate tooth profile metric for computation 
of the Fourier-series coefficients of the mean deviation trans- 
mission error components in the case of spur gears. See Fig. 4 
of Ref . 1 . 

Fourier-series coefficients and power spectra of random 
deviation components . In the case of spur gears, Eq . (55) 
expresses the Fourier-series coefficients of the random deviation 
transmission error component from gear (•) in terms of the mesh 

transfer functions <t>. „ (n/N^ * ^A) and the discrete Fourier trans- 
( . ) 0>L (•) 
forms (n) of the expansion coefficients b j as defined by 

Eq. (53). In this case, we see from Eqs . (8l) and (82) that the 

expansion coefficients b^*^„ characterize the axial average 

( • ) 0 > 0 * ( • ) / X 

e' ' (z) of the random deviation component e ' . (y,z) from tooth 
uj u j 

face j on gear (•) — i.e.. 


40 



e^(y,z)dy . 


(105) 


e cj >(z) 4 


P/2 

1 

-P/2 


That is, from Eqs. (8l) and (82) we may show in a manner com- 
pletely analogous to that used above for the mean deviation 

component that eL^(z) is given in terms of the expansion coeffi- 
( • ) 

cients bj ' £ by 


4/ (z > - J bj ( ;^(21W)' 5 V 2z/D) , 


( 106 ) 


(•) 


where the b. Q ^ are obtained using 


(•) _ (2£+l/ 

j,oA D 


, D/2 


I e cj 


(z) P £ (2z/D)dz , 


-D/2 


(107) 


( • ) 

and where the b. 0 satisfy 

J ? 0 * 


1 

D 


D/2 


-D/2 


L e cj ,(z) 


dz = l (b i o£ )2 
£=0 Jj 


( 108 ) 


** 


F ( ’ 

Cj 


(z) Of 


Hence, in the case of spur gears, the axial averages 

the random deviation components £cj^(y,z) the tooth faces on 
gear (•) are the appropriate metrics for computation of the 
Fourier-series coefficients of the rotational harmonics of the 
transmission error from gear (♦). 


Equations (137) to (140) of Ref. 1 express the discrete 

i ( * ) i ^ 

power spectrum |a rn | of the random deviation transmission 
error component for spur gears as well as for helical gears. 
However, in the case of spur gears, the one-dimensional sequence 

of expansion coefficients bj*^, £ = 0,1,2, ••• characterizes the 
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errors on tooth j as noted above; hence, for spur gears we may 
directly replace the index m in Eqs . (137) to (140) of Ref. 1 
by the index l. For a specific error component designated by 
the index a, Eqs. (84) to (87) apply to the spur gear case also 
by setting k = 0 as above. As in the helical gear case, when 
several random error components, 0£ are present, cross-spectrum 

terms must be included to obtain all contributions to | ^ I 2 
as indicated by Eq. (140) of Ref. 1 unless these various random 
error components are uncorrelated. 

Modified forms of mesh transfer functions . In the case of 
spur gears, the first form of our modified mesh transfer function 
for tooth-spacing errors is obtained from Eq. (89) by setting 
(L/A) = 0: 



except 

n*=0 


where the rect function is defined by Eq. (90). As in the case 
of Eq. (89), Eq. (109) will display small discontinuities at 

values of n in the neighborhoods of | (n/N^ ^)-n’| = ^ for all 

n*. These discontinuities will be more pronounced for spur gears 
than for helical gears. To remove these discontinuities, we 
can substitute for the rect function in Eq. (109) the weighting 
function of Eq. (91). When this substitution is made, and when 
appropriate adjustment is made for the behavior of $k£,(s) near 
g = 0 as in the helical gear case, our general form of modified 
mesh transfer function given by Eq. (95) reduces in the case of 
spur gears [k=0 and (L/A)=0] to 

$ #£ (g) ~ (-i) £ (2£+l)^j ^(TTgL) 

00 

-l (w 

n *=-» 

except 
n *=0 


o l 


( Ag,n 1 


)j 0 (n 


'irL/A) j ( Ag-n * )ttL/ A] }j 


(110a) 


42 



and 


4> kJi (g) = °> k fi 0 , (110b) 


where the weighting function w 0 £(Ag,n') is defined by Eq. (94). 
Equation (110) is the preferred form of $ 0 $,(g) applicable for 
rough calculations involving spur gears for all situations 
requiring the mesh transfer functions. 


APPLICATION 

We now illustrate the above techniques by predicting the 
various components of the static transmission error spectrum from 
a set of measurements made on a mating pair of spur gears. 

Effects of variations in the axial contact ratios of helical 
gears also are illustrated. Gear parameters required in the 
computations are listed below in Table II. 


TABLE II. - GEAR PARAMETERS 



Gear ( 1 ) 

Gear (2) 

Number of Teeth 

Pressure Angle 

Outside Diameter 

Pitch Diameter 

Base Circle 
Diameter 

50 

0.3491 rad (20°) 
0.1328 m (5.230 in.) 
0.1270 m (5.000 in.) 

0.1193 m (4.698 in.) 

53 

0.3491 rad (20°) 
0.1405 m (5.530 in. ) 
0.1346 m (5.300 in. ) 

0.1265 m (4.980 in. ) 


From the above parameters, we must compute L, A, and D. 

From Eqs . (Dll) and (D12) of Ref. 1, we find = 0.007457 m 

(0.2936 in.), L (2) = 0.007503 m (0.2954 in.); hence, 

L' = ^+L' ^ = 0.01496 m (0.5890 in.). However, the above 

computation does not account for tip rounding. Subtracting 
0.00132 m (0.0520 in.) from the above value of L' to get the 
true value of path of contact length L illustrated in Fig. D-l 
of Ref. 1, we obtain 
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L = 0.01364 m (0.5370 in.) . 


(Ill) 


The base pitch A is 27rR|- ) /N where is the base circle radius and 
N is the number of teeth — both quantities taken on the same 
gear. Using the above parameters from either gear, we find 

A = 0.007498 m (0.2952 in.). (112) 

Finally, using Eq. (C8) of Ref. 1, D = Lsin<f>, where <J> is the 
pressure angle, we find 

D = 0.004665 m (0.1837 in.). (113) 

Mesh-Attenuated Tooth-Spacing Error Spectra 

Figures 3 to 14 display the steps involved in computing the 

discrete power spectra |a^^|oo from the tooth-spacing errors of 

gears (1) and (2), where |a£j^loo is defined by Eq. ( 87 ) with 
k = 0 and i = 0 . 


Tooth-spacing errors 

Figures 3 and 6 show the tooth-spacing errors of gears (1) 
and (2). The abscissa of Figs. 3 and 6 is tooth number j and 
the ordinate is tooth-spacing error measured in micrometers 
(mu m) . Two sets of data are shown in each figure. The lower 
curve in each figure displays the accumulated tooth-spacing 
errors, which are the errors in tooth-face positions measured 
as a function of tooth number on the base circle from the posi- 
tion of an arbitrary tooth designated as tooth number zero. The 
lower curves are discontinuous because only one accumulated 
tooth-spacing error is associated with each tooth on a gear. 

Notice - e.g., in the case of Fig. 3, which applies to the 50 
tooth gear — that tooth numbers 0 and 50 have the same accumulated 
errors because they are the same tooth. The upper curves in 
Figs. 3 and 6 are tooth-to-tooth errors in tooth spacing. 
Therefore, these upper curves display the magnitudes of the 
jumps in the lower curves. The accumulated error curves shown 
in Figs. 3 and 6 are computer plots of similar curves that were 
obtained from measurements made on gears (1) and (2). Procedures 
for making such measurements are well known — e.g., see pp. 77-78 
of Ref . 6 . 
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FIGURE 3. TOOTH-SPACING ERRORS OF GEAR (1). 
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FIGURE 5. TOOTH-SPACING ERROR rms HARMONIC AMPLITUDES OF GEAR (1). 
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In Pigs. 3 and 6 , the origins of the ordinates of the accumu- 
lated error plots were chosen so that the average accumulated 
tooth-spacing error on each gear is zero. Thus, these plots 
represent random deviation components. In fact, if we interpret 
the tooth-spacing errors as measurements of the errors in the 
mean positions of the tooth faces, then we may identify the 
accumulated error of a generic tooth j on gear (•) with the 

expansion coefficient bj *] 0 - see Eqs. (82) and (107) and recall 
that Po(x) = 1. Also see Eq. (50) of Ref. 1 . Hence, the accu- 
mulated error plots in Figs. 3 and 6 satisfy 





(• 

j , 0 


0 


(114) 


Because of the circular nature of the gears, we must define the 
tooth-spacing errors so that 


b 


(•) 
p N 


(•) 


+j jO o 


(• 

j,°0> 


P = 0 , +1 , + 2 , • • • . 


(115) 


Equations (84) to (87) can be used with k = 0 , l = 0 to com- 
pute the tooth-spacing error contribution |a^^|j 0 to the discrete 
spectrum of the static transmission error of gear (•). To carry 
out this computation, we require the mean-square tooth-spacing 
errors defined by Eq. (84) with k = 0, l = 0. These values were 
computed from the accumulated errors shown in Figs. 3 and 6 and 
are 


(b 


( 1 ) 

o o 


2 


16.69 (Mu meters) 2 . 


(b^ ; ) = 19.17 (Mu meters) 2 ; 

( 116a, b ) 


hence, the corresponding rms accumulated errors are 4.085 and 
4.378 mu meters respectively, which are consistent with the 
general appearance of Figs. 3 and 6 . 

Tooth-spacing error autocorrelation functions 

We next require the normalized autocorrelation functions 
P 0 o (J ) of the tooth-spacing errors, which are defined in terms 
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of the accumulated errors bj jj q by Eq. (85) with k = 0, l = 0. 

Figures 4 and 7 show these autocorrealtlon functions for gears 
(1) and (2) respectively. The ordinate of Figs. 4 and 7 is 

and the abscissa is tooth number difference j, as may 
be seen from Eq. (85). When the tooth number difference is zero, 

we may see from Eqs. (84) and (85) that p[ 0 ^(O) = 1. Moreover, 

because tooth number ^ is the same tooth as tooth number 0, 

we must have Pg’^(N^) = 1, as we may also see from Figs. 4 and 

7. Although the abscissas in Figs. 4 and 7 span only the range 

of values 0 < j < p^(j) must satisfy 

p^tpN^+j) = p 0 ( ; } (J)> P = 0, +1, +2,..* (117) 

because of the circular nature of the gears as indicated by Eq. 
(115) • Furthermore, from Eq. (85) we see that 


P 


(•) 
o o 


(-j) 




( 118 ) 


which when combined with Eq. (117) yields the symmetric character 

of Figs. 4 and 7- From Eqs. (85) and (114), we see that Po^(j) 
also must satisfy 


N 


C^T 


* ^-1 

l 

j=0 


P ( oo )(J) 


= 0. 


(119) 


The properties of Eqs. (114), (115), and (117) to (119) apply to 
all classes of random errors designated by the dual indices k£ 
used in Eqs. (84) to (87). 

Tooth-spacing error spectra 

According to Eq. (87) applied to the case of tooth-spacing 

errors (k=0, H=0) , we require (b£* ^ ) Sg*^(n) to characterize the 
tooth-spacing errors of gears (1) and (2) for later computation 
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of the mesh-attenuated error spectrum |a^ ^| 2 . 

l a rn ; loo ls the tooth-spacing error contribution 
spectrum, as may be seen from Eqs. (109), (110), 
Ref. 1, we have chosen to compute here for gears 
one-sided tooth-spacing error spectrum 


However, since 
to a two-sided 


and (140) of 
(1) and (2) the 


G^’tn) 4 2 (b';V s' # °(n) , (120) 

where S^ n (n) is the finite discrete Fourier transform [22] of 

( . ) 

the normalized autocorrelation function p ( j ) as defined by 

Eq. (86) for k = 0 and l = 0. The square roots of G^*' (n) for 
gears (1) and (2) are plotted in Pigs. 5 and 8. Hence, the 
ordinate in Pigs. 5 and 8 is a direct measure of the rms values 
of the various harmonic components of the tooth-spacing errors. 

To provide the capability for displaying a large range of values, 

( • ) J' 

we have plotted [Gq 0 (n)] 2 on a logarithmic scale of ordinates 
labeled in units of micrometers (mu m) . The abscissa of Figs. 

5 and 8 is rotational harmonic number n. 

Figures 5 and 8 both are line spectra with fundamental har- 
monic number n = 1 corresponding to the fundamental rotational 

( • ) 

harmonic of each gear, which has a period equal to N A, the 
circumference of the base circle of gear ( • ) . The spectra shown 
in Pigs. 5 and 8 both are periodic in harmonic number with period 

equal to the number of teeth N^*^, which is = 50 in the case 

(2) 

of Fig. 5 and N = 53 in the case of Pig. 8. Two periods are 
shown in each figure. This periodic behavior is a general 

property of the spectra s£o^(n) defined by Eq . (86), which are 

^ } 
periodic in n with period equal to the number of teeth N v * ; . 

Periodic behavior of this type is a well-known property of the 
finite discrete Fourier transform [J7] and follows from the fact 
that, for integer p, we have 

exp[-i2ir (n+pN^ * ^ ) j/N^ * ^ ] = exp(-i27mj/N < ‘ * ^ ) , (121) 

which when substituted into Eq. (86) yields 

S kI^ n+pN ^*^ = S k£^ n ^ * (122) 
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From Figs. 5 and 8, we notice that the harmonics located at 
( • ) 

n = pN , p = 0, ±1, ±2,*** are absent. This absence is a con- 
sequence of the periodic behavior of s£j^(n) described above, 
and the fact that for n = 0 we have from Eq. (86), 


S M >(0) 


N 


FT 


’*-1 

l 

3=0 


><•> 

\SL 


(j ) 


(123) 


which is zero for tooth-spacing errors, as may be seen from Eq. 
(119). This result is a special case of the general result 
proved in Appendix H of Ref. 1 - i.e., that the rotational har- 
monics from all random components of the static transmission 

( • ) 

error are zero at values of n = pN , p = 0, +1, ±2,***, which 
are the fundamental and multiples of the tooth-meshing harmonics. 

From Figs. 5 and 8, we notice that the fundamental harmonic 
n = 1 is nearly an order of magnitude larger than all other har- 

. . ( . ) 

monies in the range [n| < N /2, which is the point of symmetry 
in the first period of the spectra. This large fundamental 
reflects the fact that the accumulated tooth-spacing error curves 
in Figs. 3 and 6 appear almost sinusoidal with a single period — 
which also results in the almost sinusoidal behavior of the auto- 
correlation functions shown in Figs. 4 and 7- The almost sinu- 
soidal appearance of the lower curves in Figs. 3 and 6 is typical 
of accumulated tooth-spacing error measurements — e.g., pp . 77-78 
of Ref. 6 — and is the result of an eccentrically located axis 
of gear rotation in the error measurement or gear cutting opera- 
tions . 

The symmetric and periodic behavior of the discrete tooth- 
spacing error spectra shown in Figs. 5 and 8 gives rise to a 
repetition of the large low-order harmonics in the neighborhoods 
of the tooth-meshing harmonics, which occur at harmonic numbers 

n = pN^*), p = 0, ±1, ±2,»»*. The repetition of these low-order 
harmonics in the neighborhoods of the first two tooth-meshing 
harmonics, which occur in Figs. 5 and 8 at n = 50 and 100, and 
n = 53 and 106 respectively, is easily seen in the figures. This 
repetitive character of the spectra shown in Figs. 5 and 8 is a 
fundamental consequence of the mathematical form of Eq. (86) — 
which is a finite discrete Fourier transform. This form of Eq. 
(86) is a consequence of the fact that gears are circular with 
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finite numbers of discrete teeth — and is not the result of an 
approximate sampling procedure, which is the usual motivation 
for representing physical phenomena by the finite discrete Fourier 
transform. The repetition of the low-order harmonics in the 
neighborhoods of the tooth-meshing harmonics gives rise to the 
so-called sidebands, which are observed in experimentally obtained 
gear system vibration spectra [25]. 

The magnitudes of the strongest lines in Figs. 5 and 8 can 
be checked approximately using the Fourier mate to the finite 
discrete Fourier transform [27] of Eq. (86): 

p ki )(J) = N J 0 1 s ki >(n) exp (^f) > ^ 


which gives, when evaluated at j = 0, 


p k i )(0) 


N ( *^-l 

I 

n=0 


s ki )(n) 


- 1 , 


(125) 


where Eqs. (84) and (85) were used in arriving at the second 
equality. Let us suppose that all of the "power" in Sg^(n) w 
in the fundamental sf ^(1) and its equal sideband S^*^(N^’^-1) 
From Eq . (125), we would have in this case (1) = 

^-1) = 1/2; hence, from Eq. (120), we have G^’^(l) = 

(•) 


G 


(N (,) -l) = (b n V) 


(•) 

00 ^ “ A ' “ ^ u oo 

tooth-spacing error. 


, which is the mean-square accumulated 
The values of [gJ/Q)]^ shown in Figs. 5 


and 8 are slightly less than the rms errors f ( b ^ ^ ^ ) 
and 4.378 Mu meters as predicted. 


of 4.085 


Figures 9 and 12 display the same infomration as Figs. 5 and 
8, except that Figs. 9 and 12 show the envelopes of the line 
spectra shown in Figs. 5 and 8 on a logarithmic scale on the 
abscissa. Slightly over ten cycles of the periodic spectra are 
shown in Figs. 9 and 12. The ordinate axes of Figs. 5> 8, 9 } 
and 12 all are plotted to the same logarithmic scale. Notice 
that Figs. 9 and 12 have discontinuous slopes with the discon- 
tinuities located at the positions of the rotational harmonics 
of each gear. 
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Mesh transfer functions for tooth-spacing errors 

According to Eq. ( 87 ), we next require the mesh transfer 

function $ 00 (n/N^ ^A) in order to compute |a£ n ^| 00 for gears (1) 

and (2). Figures 10 and 13, which are nominally identical, 

display the magnitudes |<j> 00 (n/N^ ^A) | of these transfer functions. 
The logarithmic scales of both axes in Figs. 10 and 13 are 
identical with those in Figs. 9 and 12. 

Equation (109) was used to compute Figs. 10 and 13. Thus, 
we observe small discontinuities in Figs. 10 and 13 approximately 
midway (on linear coordinates) between rotational harmonics that 

( • ) 

are integral multiples of N — as expected from the discussion 
of Eq . (109). These discontinuities would not have been present 

if Eq. (110) had been used to evaluate 4> 00 (n/N v 'A). 

As n tends to zero. Figs. 10 and 13 show that | (n/N^ * ^A) | 
tends to unity as required by Eq. ( 78 ). In addition, we see by 

comparing Figs. 9 and 12 with Figs. 10 and 13 respectively, that 
~ ( . ) 

<J> 00 (n/N A) is zero at the positions of the tooth-meshing 
harmonic numbers, which occur at rotational harmonic numbers 

that are integral multiples of ^ — i.e., (n/N^ ^A) = p/A 
where p = ±1, ±2,***. 

Examination of Figs. 10 and 13 shows that for large n, the 

envelope of the maxima of |<f> 00 (n/N A) | decays linearly with 
increasing n. The maxima are controlled by the first term in 
the right-hand side of Eq. (109): 


j 




sin[(n/N ( * ^ttL/A] 
(n/N^ * ^ )ttL/A 


( 126 ) 


( * ) 

Hence, for large n/N , the envelope of the maxima of 

|<j> Q0 (n/N^ ^A) | has a slope of -1 on logarithmic coordinates, 
which is equivalent to -6.02 dB per octave. In fact, the envelope 

of the maxima of |$ Q (n/N^*^A)| is very similar to the behavior 
of a first-order system — e.g., see pp. 12-17 of Ref. 19 — which 
can be approximated by two straight lines on logarithmic 
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coordinates. This approximation of the envelope of the maxima of 
|$ oo (n/N ( ^ A ) | in the case of spur gears is 


, (n/N^ * ^ ) < (irL/A) -1 


|$ 00 (n/N ( * } A) 


'env 


([ (ttL/ A)n/N^ * ^ ] 1 , (n/N^ * ^ ) > (irL/A) -1 , (127) 
where we see that (n/N^) = (ttL/A ) -1 locates the breakpoint. 


Mesh-attenuated tooth-spaaing error spectra 

When the tooth-spacing error spectrum Gg 0 ^(n) of Eq. (120) 
is multiplied by the squared mesh transfer function 

|$ (n/N^ } A) | 2 , we obtain the one-sided mesh-attenuated tooth- 

spacing error spectrum 2|a^ n | 00 — according to Eqs. ( 87 ) and 
(120). The envelopes of the square roots of these one-sided 
line spectra are plotted in Figs. 11 and 14 for gears (1) and (2) 
respectively. A line occurs at each discontinuity in slope in 
the two figures. Figures 11 and 14 are, respectively, the 
products of the functions plotted in Figs. 9 and 10 and Figs. 

12 and 13 . The same logarithmic scales are used in all of these 
plots . 

Since the mesh transfer functions in Figs. 10 and 13 are 
essentially unity at the positions of the fundamental rotational 
harmonic, the fundamental rotational harmonics in Figs. 9 and 
11 have the same value as do those in Figs. 12 and lh. Thus, 
the meshing action of the gears provides no attenuation of the 
fundamental rotational harmonic of tooth-spacing errors. However, 
the higher-order harmonics receive considerable attenuation — as 
may be seen by comparing Figs. 9 and 11 or Figs. 12 and 14. The 
asymptotic slopes of the envelopes of the maxima in Figs. 11 and 
14 is -1, which corresponds to -6.02 dB per octave — as was the 
case for the mesh transfer functions. 

It is particularly interesting to compare the strengths of 
the sidebands nearest to the tooth-meshing harmonics in Figs. 9 
and 12 with the same mesh-attenuated sidebands appearing in 
Figs. 11 and 14. These mesh-attenuated sidebands are strongly 
attenuated in comparison with their neighboring harmonics, which 
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show much less attenuation. This strong attenuation occurs 
because the sidebands closest to the tooth-meshing harmonics 
partially fall into the nulls of the mesh transfer functions 
that occur at the tooth-meshing harmonics, as shown in Pigs. 10 
and 13 . 


Tooth-Meshing Harmonic Contributions and Effects of 
Different Axial Contact Ratios 

Figures 15 to 17 illustrate contributions to the tooth- 
meshing harmonics of the static transmission error and Pigs. 18 
to 20 illustrate effects of different axial contact ratios for 
the case of tooth-spacing errors on helical gears. 

Load-de-pendent component 

In the case of spur gears with constant tooth-pair stiffnesses 
per unit length of line of contact, Eq. (98b) gives the Fourier- 
series coefficients of the tooth deformation component of the 
static transmission error for a pair of meshing gears. Figure 
15 shows the normalized values of the corresponding rms ampli- 
tudes of these tooth deformation harmonics: 


✓2 


a Wn 

— Jo 

sin(mrL/A) 

w„/k t 

— Yd 

nuL/A 


( 128 ) 


Comparing the positions of the harmonics in Fig. 15 with Figs. 11 
and 14, we see that the tooth-meshing harmonics of Fig. 15 occur 
at the positions of the nulls shown in Figs. 11 and 14, as noted 
earlier. Since the magnitudes of the harmonics shown in Fig. 15 
are normalized by the mean deformation W 0 /Krp, no direct compari- 
sons should be made between the heights of the lines in Fig. 15 
and the envelopes shown in Figs. 11 and 14. 


Mean deviation components 

Equation (50) is our general expression for the Fourier- 
series coefficients of the mean deviation component of the trans- 
mission error from gear (•). In the present case, where we have 
assumed the tooth-pair stiffness to be constant, the expansion 

coefficients a.o^ for use in Eq. (50) are given by Eq. (102), 

/ 9 \ U 

where m^* '(z) is the axial average of the mean deviation 

mc*)(y,z) of the tooth-faces of gear (•) from perfect involute 
surfaces — as may be seen from Eq. (100). 


64 



MESH -ATTENUATED rms AMPLITUDE 












MESH -ATTENUATED rms AMPLITUDE 





sj 



FIGURE 18. TOOTH-SPACING ERR' 
WITH (L/A) = 1.81! 
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FIGURE 19. TOOTH-SPACING ERROR MESH TRANSFER FUNCTION OF HELICAL GEAR 
WITH (L/A) = 1.819 AND (FL/AA) = 3.19. 
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FIGURE 20. TOOTH-SPACING ERROR MESH TRANSFER FUNCTIO'! OF HELICAL GEAR 
WITH (L/A) = 1.819 AND (FL/AA) = 9.251. 
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Pour measurements of the deviations of the tooth profiles 
of the 50 tooth gear from perfect involute profiles are shown 
in Fig. 21. The four teeth measured were picked at approximately 
90° intervals around the gear. The tips of the teeth are to the 
right and the roots are to the left. Distances along the 
abscissa in Pig. 21 are measured in the same units as the 
coordinate s in Ref. 1. Thus, the region of tooth contact along 
the abscissa in Pig. 21 is of length equal to the path of contact 
length given by Eq. (Ill) - as indicated in Pig. 21. Notice that 
one-half of the value cited earlier for tip rounding is accounted 
for in determining the region of contact in Pig. 21. The other 
one-half is accounted for in determining the region of contact 
on the teeth of the mating gear. 

The variable of integration z in Eq. (102) is related to 
coordinate s by Eqs . (16) and (17) of Ref. 1, where y = 0 for 
spur gears since ip = 0 in this case - see Eq. (18) of Ref. 1. 

Thus, for spur gears, we have z = s sin<J>; hence, D = Lsin<f>, as 
may be seen from Eq. (C8) of Ref. 1 where D is the active tooth 
depth illustrated in Pig. 4 of Ref. 1: D = 0.004665 m (0.1837 

in.). The abscissa in Pig. 21 therefore is proportional to the 
variable of integration z with origin centered in the region of 
contact, which is of width D when measured in units of z. 

Values of the ordinates in Fig. 21 are deviations of the 
tooth profiles from perfect involute curves as defined in Fig. 3 
of Ref. 1 and the accompanying discussion. Since these deviations 
are defined as positive when they are "equivalent" to removal of 
material from perfect involute tooth faces, the positive direction 
of the ordinate in Fig. 21 is downward. The location of the 
origin of the ordinate in Pig. 21 is irrelevant since the origin 
has an effect only on the zeroth order expansion coefficients 
( • ) 

a 00 determined by Eq. (102), and we have seen that the corre- 
sponding expansion function $ (n/A) in Eq. (50) is zero except 

for n = 0 — see Eq. (78). JJ 

If we interpret each of the four profile deviation measure- 
ments in Fig. 21 as a measurement of the axial average of the 
tooth-face deviation, then the average value. of these four 
curves (after removal of vertical separations) becomes our best 

estimate of the mean profile deviation mQ^(z) for the 50-tooth 
gear. However, all four curves shown in Fig. 21 are virtually 
identical except for their vertical translations; hence, we 

computed the expansion coefficients a^ given by Eq. (102) from 
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FIGURE 21. MEASURED DEVIATIONS FROM INVOLUTE PROFILES OF 4 TEETH 
SPACED APPROXIMATELY 90° APART ON GEAR (1). 



only one curve — the second curve from the top shown in Pig. 21. 
This computation was begun by first reading 54 equi-spaced points 
off the second curve. To accommodate the fluctuations of the 
P£(2 z/D) in Eq. (102), we then used linear interpolation to obtain 
five additional points on the second curve between each of the 
original 54. The integrals in Eq. (102) were then evaluated 
using the trapezoidal rule operating on the integrand evaluated 

at the resulting 319 points. Expansion coefficients aj;]^ were 
computed through £ = 20 in this manner. As a check, the profile 
was then reconstructed from these expansion coefficients using 
Eq. (101) — a fit to the original profile with negligible error 
was obtained. The expansion coefficients through £ = 20 were 
then used together with Eq. (99a) to calculate the Fourier-series 

coefficients using Eq. (50). The corresponding rms ampli- 

tudes., which are /2|a^^|j are plotted in Pig. 16. 

The same procedure was used with the mating gear (2), which 
has 53 teeth. The only difference in the procedure for this 
gear was our reversal of the direction of the z-axis — i.e., 
the tooth tips were on the negative side of the origin of the 
abscissa — because the tips of gear (2) mesh with the roots of 

gear (1). The rms amplitudes, which are /2’|a I f in ; | for gear (2), 
are shown in Pig. 17. The harmonics shown in Figs. 15, 16, and 
17 all occur at the same positions, which are positions of nulls 
in Pigs . 11 and 14 . 

The total contribution of the mean profile deviation com- 
ponents from gears (1) and (2) would be obtained by adding the 

complex Fourier-series coefficients and from gears 

(1) and (2). Finally, the Fourier-series coefficients of all 
components of the tooth-meshing harmonics would be obtained by 
adding to the coefficients of the mean profile deviation com- 
ponents the complex Fourier-series coefficients of the load- 
dependent component as indicated by Eq. (57) or Eq. (65) of 
Ref. 1. 

Mesh transfer functions for tooth-spacing errors on helical gears 

So far, we have illustrated our results for spur gears only. 
Let us now compare some of these results with results that would 
be predicted for helical gears. We recall that the tooth-spacing 
error spectra shown in Figs. 5, 8, 9, and 12 are independent of 
all gear design parameters except number of teeth. Hence, these 
spectra could just as easily have been computed from measurements 
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made on helical gears. It follows that we can estimate differ- 
ences in the mesh-attenuated tooth-spacing error spectra that 
would be observed in comparing spur and helical gears — and 
differences that would be observed in comparing helical gears 
with different axial contact ratios — by comparing differences 
in the mesh transfer functions for tooth-spacing errors between 
spur and helical gears, and among helical gears with different 
axial contact ratios. 

Figures 18 to 20 show the magnitude of the tooth-spacing 

error mesh transfer functions I $ 0 0 (n/N^ ^A) | for 50 -tooth helical 
gears with three different values (1.1, 3.19, and 9.251) of the 
axial contact ratio FL/AA. Each of Figs. 18, 19, and 20 was 
computed by Eq. (88) using k = 0, SL = 0 and the same value of 
transverse contact ratio, (L/A) = 1.819, that was used to compute 
the spur gear mesh transfer functions shown in Figs. 10 and 13. 
Furthermore, the logarithmic scales on the ordinates and abscissas 
of Figs. 18 to 20 are identical with those used in Fig. 10. Thus, 
Figs. 10, 18, 19, and 20 provide a direct illustration of the 
very substantial decrease in the tooth-spacing error vibratory 
excitation that is provided by helical gears in comparison with 
spur gears. The attenuation of the excitation is increased with 
increasing axial contact ratio, and the resulting improvements 
increase with increasing harmonic number. Notice the very sub- 
stantial increase in attenuation in the higher harmonics predicted 
for the helical gear with the relatively small axial contact ratio 
of (FL/AA) =1.1 shown in Fig. 18 in comparison with the spur 
gear result shown in Fig. 10. Finally, we see that the mesh 
transfer functions shown in Figs. 18 to 20 all have nulls at the 

tooth meshing harmonics of (n/N^^A) = p/A for p = ±1, ±2,***. 

The tickmarks on the abscissas in the figures show the locations 
of the first ten tooth-meshing harmonics. 

We may approximate the behavior of the envelopes of the 
maxima in Figs. 18 to 20 from the first term in Eq. (88) for 
k - 0, l =0: 

/ n ttFL\ . / n ttL\ = sin[(n/N ( * ] )ttFL/AA] sin[ (n/N ( ‘ \ )irL/A ] 

°\N * ^ AA / J °\N ( * ) A / (n/N ( ’^ttFL/AA (n/N ( * } )ttL/A 


(129) 

Hence, for large n/N^*^, the envelope of the maxima of 
| $ 0 (n/N^ ’ ^ A) | has a slope of -2 on logarithmic coordinates. 
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which is equivalent to -12.04 dB per octave. This behavior is 
typical of second-order systems - e.g., see pp. 24-29 of Ref. 

19 - and can be approximated by two straight lines on logarithmic 
coordinates. The resulting approximation of the envelope of the 

maxima of |<f> 00 (n/N^ ^A)| -in the ease of hetioat gears is 



where we see that (n/N^) = [it 2 (FL/AA ) (L/A ) locates the 

breakpoint . 

Mesh Transfer Functions for Higher-Order Error Contributions 

Figures 18 to 20 have illustrated the effects of different 
axial contact ratios of helical gears for one particular error 
component — the tooth-spacing component k = 0, Z = 0. Figures 
22 to 27 illustrate the behavior of the magnitudes of the mesh 
transfer functions of different error components while holding 
the axial and transverse contact ratios at constant values: 


FL 

AA 


3.19, 


1= 1.819 


(13D 


The order of the indices k,£ in going from Figs. 22 to 27 is that 
indicated by Fig. 1. The logarithmic scales of the ordinates and 
abscissas in Figs. 22 to 27 are identical with those used in 
Figs. 18 to 20. 

Equation (88) was used to compute the values of 

l$k*,(n/N^A)| shown in Figs. 22 to 27. Thus, Fig. 22 is a 

repeat of Fig. 19 - and shows the nulls of $ 00 (n/N^A) that 
occur at the tooth-meshing harmonics as discussed earlier. The 
other mesh transfer functions shown in Figs. 23 to 27 do not 
possess these nulls. However, the envelopes of the peaks of 
all of the mesh transfer functions have the same asymptotic 
slope (-12.04 dB per octave) in the large harmonic number region. 
Their behavior differs in the low harmonic number region. 
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MESH TRANSFER FUNCTION, <f> nn ( n /N l ' 'A ) 








MESH TRANSFER FUNCTION, <f> 0] (n/N m A) 



ROTATIONAL HARMONIC NUMBER (n) 


FIGURE 23. INVOLUTE SLOPE DEVIATION (k=0, 1 = 1 ) MESH TRANSFER FUNCTION 
OF HELICAL GEAR WITH (L/A) = 1.819 AND (FL/AA) = 3,19. 



MESH TRANSFER FUNCTION, |$, 0 (n/N ln A) 
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FIGURE 24. LEAD MISMATCH DEVIATION (k=l, *=Q) MESH TRANSFER FUNC 
OF HELICAL GEAR WITH (L/A) = 1.819. AND (FL/AA) = 3.19 
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Envelopes of mesh transfer functions for helical gears 


We may quantify the envelope behavior of the maxima in Pigs. 
22 to 27 using the magnitude of the first term in Eq. (88): 

Prom Eq. (36) and the above expression, we see that the envelope 
of the maxima in Figs. 22 to 27 must asymptotically approach for 

large n/N ^ the form 


l$ k 5, (n/N( * )A) 


env 


[ ( 2k+l ) ( 2&+1 ) ]^ / n V 
tt 2 (FL/AA) (L/A) \N^‘ V 


( 132 ) 


which is independent of k and l except for the term in the numera- 
tor. For small x, the first term in the Maclaurin series expan- 
sion of J m (x) gives - see, e.g., p. 437 of Ref. 4 - 


m 


j„( x ) = 

m 


x 


1*3*5' 


( 2m+l ) 


m+2. 

0(x ) 


(133) 


which we may use to evaluate the magnitude of the first term in 
Eq. (88) to yield the behavior of |$ k ^(n/N^ ^)| env /or small 
n/ »<'>: 


l* k , (n/N ( ’ ' a) | 

KX ' env 


TT k+g, [(2k+l) (2&+l)] ig 

[ 1 • 3 * 5 * * * ( 2k+l ) ] [ 1 • 3 * 5 * * * ( 2 A+l ) ] 



(134) 


Hence, for small n/N^ \ ^^^(n/N^ ^A) 
coordinates with slope equal to k+£. 

zero slope for small n/N^*^ observed in Fig. 22; it also predicts 
that Figs. 23 and 24 should have the same asymptotic slope for small 


is linear on logarithmic 
This behavior predicts the 


n/N^*\ and that Figs. 25, 26, and 27 all should have the same 

somewhat larger value of slope for small n/N^ ^ — which is con- 
sistent with the behavior shown in the figures. 
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We see from Pigs. 22 to 27 that we can obtain a reasonably 

good approximation to the envelope of the maxima of |<j>. 0 (n/N^A)| 

for all values of n/IT ' by extending the large and small n/IT * ; 
approximations given by Eqs. (132) and (13*0 respectively, until 

they meet. The breakpoint (n/N^*^) bp is obtained by equating 

the right-hand sides of Eqs. (132) and (134) and solving for 

n/IT 

1 

[1*3 *5* • • (2k+l ) ][1 *3 * 5 * » • (2ft+l ) ] | 2+k+£ 

(PL/AA) 1+k (L/A) 1+Jl ) * ^35) 


Prom Eqs. (75) and (83), we see that the squares of the 

( • ) ( • ) 

expansion coefficients a£j/ and bj , are direct measures of the 
mean-square contributions of the virious terms k£ in our expan- 
sions of the mean and random tooth-face deviations mc*^(y,z) and 
e^ } (y,z). Hence, from Eq. (48) and Eqs. (84) to (87), it is 

apparent that the mesh transfer functions may be compared 

with one another to determine the relative attenuation provided 
by the meshing action of a gear pair on the various classes of 
errors defined by the dual index k£. In the large harmonic 
number region, we see from Eq. (132) that - except for tooth- 
spacing errors in the immediate vicinity of the tooth-meshing 
harmonics as noted earlier — the only systematic dependence of 

| $. 0 (n/N ^A) I on k and % is the factor [ ( 2k+l ) ( 2&+1 ) in 

k£ 'env 

the numerator of Eq. (132). Hence, the mesh provides somewhat 
less attenuation for the higher-order terms k£ than for the 
low-order terms. The differences in the low harmonic number 
region are more striking as k and % are varied. In particular, 

we see from Eq. (134) that at (n/N^) = 0, all of the mesh 
transfer functions are zero except that associated with tooth- 
spacing errors, k = 0, l = 0. 

Envelopes of mesh transfer functions for spur gears 

To develop expressions for the envelopes of the mesh transfer 
functions for spur gears that are comparable to Eqs. (132), (134), 
and (135) — which apply to helical gears — we require the magni- 
tude of the first term in Eq. (99a) evaluated at g = n/N^A: 



1 

IT 
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(24+1)' 


/ n irL \ 

It 77 4 ) 


Using Eq. (36), we see from the above expression that the envelope 

of the maxima of | $ 0 (n/N^ * ^A) | must asymptotically approach for 
( .) * 
large n/N the form 


l$ 0 <,(n/N (-) A)| 

env 


{21+1)' 2 

tt(l/A) 



(136) 


which corresponds to a slope of -6.02 dB per octave. We can use Eq. 
(133) to yield the behavior of I <j> ft 0 (n/N^ * ^ A) I for small n/N (,) 

OX/ 0X1 V 


|$ oA (n/N (,) A) 


env 


tt & (2£+1)^ 

[ 1 * 3 * 5 * * * ( 2SL+1 ) J 



(137) 


The breakpoint for spur gears (n/N^) b is obtained by equating the 
right-hand sides, of Eqs. (136) and (137) and solving for n/N^*\ 



1 

[1*3*5* * *(2&+l)] 1+Jl 
tt (L/A ) 


( 138 ) 


which is used to divide the regions of application of Eqs. (136) 
and (137). 


We may solve for the maximum value of the envelope 
|$ oA (n/N^ * ^A) | env max by evaluating Eq. (136) or Eq. (137) at 

(n/N (<) ) = (n/N (,) ) bp : 


|$ o£ (n/V * \) 


env max 


[1*3*5 


( 2 &+ 1 ) 3 * 

• • (24+1)] 


1 / ( 1 + 4 ) 


(139) 


which is independent of L/A. Table III shows the dependence 
of (n/N^*^) bp on H computed from Eq. (138) and lists values 
|$ o £ (n/N^A) | env max computed from Eq. (139). We 
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TABLE III. - ENVELOPE BREAKPOINT POSITIONS AND MAXIMA FOR 
SPUR GEAR MESH TRANSFER FUNCTIONS 


l 

0 

1 

2 

4 

8 

16 

32 

ttL / n \ 

1 

1.73 

2.47 

3.94 

6.88 

12.76 

24.54 


X 

I $ . 0 ( n/N ^ ^A) 

1 'env max 

1 

1 

0.907 

0.762 

0.599 

0.45 

0.329 


see from the table that for any fixed value of L/A the break- 

( • ) 

point value of n/N increases monotonically with l while the 
value of the envelope maximum decreases monotonically with l . 

Since the magnitude of the first term in Eq. (88) given by 
the expression preceding- Eq. (132) has the form of the product 
of two terms each having the form of the expression preceding 
Eq. (136), it follows from the behavior shown in Table 3 that 
in the case of helical gears the maximum value of the envelope of 

l^g^n/N^ ^A) | will generally decrease monotonically in both k 

and SL and is strictly bounded above by the values shown in the 
table where SL should be replaced in this interpretation by the 
larger value of k or l. 
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APPENDIX A: APPROXIMATE TRUNCATION ERRORS IN EXPANSION OF 

RECIPROCAL MESH STIFFNESS 


In Appendix F of Ref. 1, an expression, Eq. (P12), was 
derived for the approximation KipjyjCx) of the reciprocal mesh 
stiffness K^ : (x) defined in Eq. (F3) of that reference: 


K 


K t [1+6K t (x)/K t ] 


(Al) 


Prom Eqs. (PI) to (P12) of Ref. 1, it follows that the approxi- 
mation K~jJj(x) can be expressed as 

/_ \ M 

K™ ( x ) = K ) l [-6 K T (x)/K ] m . ( A2 ) 

V / m=0 

The metric of error that we shall use is the fractional error 

in our estimate KipjJj(x) of K^ 1 (x). According to Eqs. (Al) and 
(A2), this fractional error can be expressed as 


K T ~ k tm = 1 


k; 


- 1 


|[1+(6K t /K t )] 



(-6K t /K t ) 



(A3) 


which is a function only of 6K T /Km and M. To a first approxi- 
mation, K t (x) is proportional to the number of teeth in contact, 
where the proportionality constant is the stiffness of one 
tooth. Furthermore, to a first approximation, we may assume 
that the number of teeth in contact fluctuates by one tooth — 
i.e., fluctuates by plus and minus one-half about the mean 
number of teeth in contact. Hence, we may approximate the 
magnitude of 6K^,/K^ by 



K,j contact ratio 


where, here, the contact ratio is interpreted as the mean number 
of teeth in contact. Table 1 was computed as a function of M 
and the contact ratio by combining Eqs. (A3) and (A4). 
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APPENDIX B: DOMINANT TERM IN SUMMATION IN EQ. (77) 

Since the numerator of each term in the summation in Eq. 
(77) is a sine function, the magnitude of the overall numerator 
for a given value of n is very sensitive to the values of L/A 
and PL/AA - but always must fall between zero and unity. Thus, 
to determine the most important term in the summation for 
general values, of L/A and PL/AA, we shall consider only the 
denominator D(n') of a typical term - i.e., 

D(n') = ( it PL/ AA ) 2 ( ttL/A ) 2 [ n ' ( n-n ’ ) 1 2 • (Bl) 

Differentiating Eq. (Bl) with respect to n' and setting the 
resulting expression equal to zero yields three stationary 
points of D(n’), which fall at values of n' = 0, n/2, and n. 
Evaluating the second derivative of D(n') at these three values 
of n* shows that D(0) and D(n) are minima, and that D(n/2) is 
a maximum point. Therefore, maximum values of [D(n')] -1 occur 
at n' = 0 and n* = n. However, from Eq. (77) we see that the 
term n* = 0 is excluded from the summation; hence, for a general 
range of values of L/A and PL/AA, we conclude that the dominant 
term in the summation is the term n* = n. 
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